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Abstract 



A comprehensive theory of nuclear magnetic relaxation in 5* = 1 Haldane gap ma- 
terials is developed using nonlinear-cr, boson and fermion models. We find that at tem- 
peratures much smaller than the lowest gap the dominant contribution to the relaxation 
rate comes from two magnon processes with T^^ ~ e~^"^^'^ , where is the smallest 
gap corresponding to a polarization direction perpendicular to the field direction. As 
the gap closes, we find that the dominant contribution comes from one magnon pro- 
cesses, and the result depends on the symmetry of the Hamiltonian. Overall the models 
agree qualitatively, except near the critical regime, where the fermion model is shown 
to be the best description. We include a thorough discussion of the effects of interchain 
couplings, nearest neighbour hyperfine interactions and crystal structure, and introduce 
a new theory of impurities corresponding to broken chain ends weakly coupled to bulk 
magnons. The work is then applied to recent measurements on NENP. We find overall 
fair agreement between available ^ data and our calculations. We finish by suggesting 
further experimental tests of our conclusions. 



ii 



Table of Contents 



Abstract 

Table of Contents 
List of Figures 

1 Introduction and Background 

[Ll] Introduction 

|1.2| Spin-wave Theory 

[T^ Non-Linear a (NLcr) Model . . . 

[T^ Boson Model 

|1.5| Fermion Model 

|1.6| Nuclear Magnetic Relaxation Rate 



2 Details of the Models 



271] NLcr Model: Temperature and Field Dependence of the Spectrum; Exact 

Results 

2.1.1| Temperature Dependence of the Gap 



2.1.21 Field Dependence of the Gap 



2.1.3 Exact Results 



|27^ The Free Boson and Fermion Models 



2.2. 1| Diagonalization 



m 



2.2.2| Discussion: Comparison of Spectra and Spin Operator Matrix El- 
ements 



3 Model Predictions for Tf ^ 

O Tf ^ foT h<^hc 



3.1.1 Isotropic Symmetry 



3.1.2| Axial Symmetry 



3.1.31 ^2 X ^2 X Z2 Symmetry 
U Close to the Critical Field . . . 
Above the Critical Field . . . . 



Summary 

4 Material Properties and Possible Effects on Experiment 

|4.1| Hyperfine Tensor 

Impurities 

Interchain Couplings 

Crystal Structure 



5 NENP: Direct Comparison with Experiment 

|5.1j The Structure of NENP and Experimental Ramifications 
|5.2| Analysis of the Data 



6 Suggested Experiments and Curious Predictions 
|6.1| Experimentally Testable Conflicts Between Models 

|6.2| Measuring Small Anisotropies 

O ESR for NENP 



7 Concluding Remarks 



IV 



Bibliography |120 



V 



List of Figures 



2.1 Free boson and free fermion dispersions witli tlie gap parameters of NENP. 

— * — * 

Top grapli: H\\b; bottom grapli: H\\c 

3.1 First non-vanisliing contribution to relaxation due to the staggered part 
of the spin 

3.2 Inter- vs Intrabranch transitions 

3.3 The gap structure for 0(3) symmetry 

3.4 The gap structure for U{1) symmetry, (a) A3 > A_l; (b) A3 < A_l . . . . 

4.1 Impurity level diagram when D' = 

4.2 Impurity level diagrams ioi D' ^ 

5.1 NENP 

5.2 Local and crystallographic axes projected onto the ac-plane in NENP . . 

5.3 A projection of the NENP unit cell onto the ac-plane showing two chains 
per unit cell 

5.4 Boson F{h, T) for fields along the b and c chain directions 

5.5 Fermion F{h,T) for fields along the b and c chain directions 

5.6 Boson F{h, T) for fields along the b and a chain directions 

5.7 Fermion F{h,T) for fields along the b and a chain directions 

5.8 Theoretical (lines) vs. experimental data (circles and squares) 

5.9 Theoretical (lines) vs. experimental data (circles and squares) 

5.10 Relaxation rate for field along the 6-axis 



vi 



5.11 r(T) — the width of the impurity resonance |107 



6.1 Dispersions for the two chain conformations and sample resonant transi- 



tions for a uniform field placed 60° from the crystallographic c-axis. . . . 116 



6.2 Resonant field vs. field orientation in the ac-plane for .19 meV transitions. |117 



vii 



Chapter 1 



Introduction and Background 

1.1 Introduction 

In 1983, Haldane derived his famous result stating that integer spin one dimensional 
Heisenberg antiferromagnets featured a gap in their low energy excitation spectrum [|l|. 
Since then, much effort has been devoted to further exploration of such systems, both 
experimentally and theoretically. The purpose of this work is to develop a theoretical 
framework for the understanding of low energy experiments on one dimensional Haldane 
gap materials. In particular, we focus on the nuclear magnetic relaxation rate, Tf^, 
although the work has relevance to many other techniques. By studying this thesis, it 
is hoped that the reader can become familiar with the tools used to understand integer 
spin Heisenberg antiferromagnetic chains with anisotropics, and can apply these tools to 
the analysis of real systems. 

There are, essentially, three important models that have so far been used to describe 
the system. In the later sections of this chapter we review the competing descriptions 
of S* = 1 antiferromagnetic spin-chains, paying some attention to their advantages and 
shortcomings. We start by outlining the traditional spin-wave theory^ used to model 
antiferromagnetism in higher dimensions. After illustrating the deficiencies in this ap- 
proach, we describe the Nonlinear a (NLa) model in some detail. This is followed by 
an analysis of a simplified yet closely related Landau- Ginsburg boson model. Last, we 
discuss a free fermion model used recently to successfully treat the case of anisotropic 
^see for a comprehensive discussfon of this topic 

1 
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spin-chains. We end the chapter with background on the nuclear magnetic relaxation 
rate, 1/Ti, for nuclear spins coupled to the spin-chain through hyperfine interactions. 

Chapter Two focuses on the details of the models, building the tools necessary for a 
detailed analysis. We discuss the temperature and magnetic field dependence of the NLcr 
model and its possible relevance to the spectrum of the boson model, as well as cite some 
exact results available in cases of high symmetry. We also diagonahze the free boson and 
fermion field theories, including on-site anisotropy effects. We derive matrix elements of 
the uniform spin operator between one particle states of magnetic excitations. These are 
used to compare the different models. 

Chapter Three explicitly describes the calculation of NMR Tf" ^, considering various 
symmetries of the Hamiltonian. We identify the leading mechanisms for low temperature 
relaxation in the presence of a magnetic field. We discuss three regimes corresponding 
to different magnitudes of the applied external magnetic field, giving expressions for the 
rate in each case. We discover that at temperatures much lower than the smallest gap 
the uniform part of the spin operator contributes most to the relaxation rate; in the 
absence of interactions, this corresponds to two magnon processes. The rate is found 
to be ~ e~^"*/^, where is the smallest gap corresponding to a polarization 
direction perpendicular to the magnetic field. As the externally applied magnetic field 
approaches a critical value, one of the gaps closes and we find the dominant process to be 
one magnon, corresponding to contributions from the staggered part of the spin. In this 
regime, we show that the fermion model is the best description and that the expression 
for Tj~^ depends on the symmetry of the Hamiltonian. 

Chapter Four deals with intrinsic effects which must be taken into account when 
analyzing experimental data. We discuss nearest neighbour hyperfine interactions; we 
show that these will contribute to order Ann/ A, the ratio of the nearest neighbour coupling 
to the local coupling. We also consider interchain couplings and show that they introduce 
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a natural infrared cutoff to the diverging density of states at the gap; for sufficiently long 
chains, they also densely fill the energy intervals between states along a finite chain. 
Finally, we introduce a new impurity theory to explain the effects of nearly free spin-| 
chain end degrees of freedom. We find that the states formed by such end spins in the 
gap, can give rise to non-trivial relaxation when coupled to the bulk excitations. 

Chapter Five applies the theory to recent experiments on the well studied material, 
Ni{C2HsN2)2N02{Cl04) (NENP). We take a close look at the crystal structure of NENP 
and identify possible terms which may be present in the Hamiltonian. We also note the 
fact, hitherto neglected, that NENP possesses two inequivalent chains in each unit cell. 
The results of Chapters Three and Four are then used to analyze experimental data. We 
find reasonable agreement for a magnetic field placed along the crystal c-axis of NENP, 
and an unexpected discrepancy for a magnetic field placed along the chain axis. The 
impurity theory is used to model low field data with qustionable results. 

The final chapter proposes further experimental tests of the theoretical predictions of 
this work. We suggest elastic neutron, electron spin resonance and further NMR studies 
to verify our own. 

1.2 Spin-wave Theory 

The Heisenberg Hamiltonian describing the isotropic antiferromagnetic spin-chain is 

H = JY^ Si - J>0 (1.1) 

i 

This arises naturally from the Hubbard model for insulators at 'half filling |^. To 
understand where this might come from, we follow the case where there is a triplet of 
possible spin states per site. On each site there are a number of valence electrons (eight 
valence electrons in the d shell of Ni^"^, for example); the degenerate electronic levels 
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are split in a way determined by Hund's rules^ and the symmetry of the crystal fields 
surrounding the ion. In some special cases (as with Ni^"^ in a field with octahedral 
symmetry), a degenerate triplet of states lies lowest. The ensuing low energy physics can 
be essentially described using effective spin 1 operators Q]. By 'half filling, we mean that 
there is an effective S = 1 triplet of states for every site in the chain (ie. there are two 
singly occupied orbitals on each site. Other orbitals are either empty or doubly occupied. 
Spins in singly occupied orbitals are aligned by Hund's rules.) Antiferro magnetism comes 
from allowing a small amplitude for nearest neighbour hopping which is highly suppressed 
by coulomb repulsion from the electrons already occupying the site. 

In the quantum case, the spin operators have the commutation relations: 

[S-, S'^] = t6,, e'^'^^r, S,- 5. = s{s + 1) (1.2) 

c 

where Sij is the Kronecker Delta Function and e"'"^ is the completely antisymmetric Levi- 
Civita symbol. 

It is easy to see that the classical Neel ground state with alternating spins is not the 
quantum ground state. To this end we write the Hamiltonian in terms of raising and 
lowering spin operators: 

= s"" ± isy 



;i.3) 



The Neel ground state is composed of spins alternating in quantum numbers between 
sites. 

|Neel >= \sl = +1, = -1, = +1, = -1 > (1.4) 



^Hund's rules maximize the total electronic spin and the total angular momentum of the electrons in 
the valence shell. 
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This state is clearly not an eigenstate of the above Hamiltonian since upon acting on 
it, the S^S~^i terms in the Hamiltonian generate states with = 0. To proceed in 
understanding the low energy properties one usually assumes that the ground state is 
approximately Neel with quantum fluctuations. The picture is that of zero point motion 
about the positions of the classical Neel spins. What we will shortly see is that the 
assumption of small fluctuations breaks down in one dimension. 

The conventional approach makes use of the Holstein-Primakov transformation. One 
begins by dividing the chain into two sublattices, "A" and "B", with adjacent sites on 
separate sublattices. On sublattice "A" one defines 

SI = s — ajai, Sj^ = a\\l2s — a\ai (1.5) 

On sublattice "B" we have 

SU, = -s + h\h, St,, = b\^2s-h\h (1.6) 

a and h are usual bosonic operators with commutation relations: 

[a, a^] = [b, = 1, [a, a] = [b, b] = (1.7) 

It can be checked that (|1.5|) and (|1.6| ) preserve the correct spin commutation relations 

— * — * 

and the constraint Si ■ Si = s{s + 1). The Neel ground state is one without bosons. 
So far no approximations have entered into the picture. However, to make progress, we 
assume that s is large. This is equivalent to a semi-classical approximation since for 
s ^ oo the commutator of the spins will have much smaller eigenvalues than the square 
of the spin variables 

[5^ S'']=tJ2 e"''-^' = 0(s) < 0{s^) (1.8) 

c 

We expand the spin operators to give 

= alVTs, Sti = b\^prs (1.9) 
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To leading order, the Hamiltonian reduces to 

H = JJ2 [-s^ + s{2alai + 2blk + aA + h^iai + h\a\ + affeti) + 0(1)] (1.10) 

i 

Fourier transforming: 



N_ 

1 



a 



e'2-^"/^afc (1.11) 



with a similar expression for 6; is the number of sites on each sublattice. Ignoring the 
constant term, we rewrite ( |1.10| ) as 



H = 2JsY^ [aW + bih + (1 + e'^'=")(afc6_, + felalj] (1.12) 



k 

N Afl 



2a is the sublattice spacing and k = nn/Na for n G [— y]. We now make the Bogoli- 
ubov transformation, 

dk = Ukb-k - Vkal (1.13) 

where, 

j^ika/2 

Uk = — 1=- (1 + csc(A;a)) ' 

vu = (-1 + csc(H)'^' (1-14) 

The d's and c's are spin wave operators corresponding to magnetic excitations (magnons) 
with = ±1 respectively. This transformation preserves the commutation relations and 
the Hamiltonian can now be written as 

H = 2Js^sin(|fca|) (4cfc + 4^) (1-15) 

k 

We see that this low energy description implies the spins are in some coherent state 
of a's and 6's built on the Neel state, but there are long wavelength Goldstone modes 
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with dispersion relation uj = 2Js\ka\ which allow each of the sublattice magnetization 
vectors to make long wavelength rotations. To ensure consistency, we now look at the 
expectation value of the magnetization (say, on the "A" sublattice), hoping to see that we 
get the semi-classical result < S'f >= s — 0(1). We invert the Bogoliubov transformation 
to get 

ttk = ulck + Vkd\ (1-16) 

and compute 

< S'- >= s- < ajai >= s - < 4«A; > 

k 

= s — a [ — |t>fcP = s — a / — (csc(A;a) — 1) (1-17) 
J 2n J An 



The last line follows from ( |1.16| ) and the fact that the true ground state has zero spin 



wave occupation number. The problem is now apparent: low wavelength modes cause 
< s — S*^ > to diverge logarithmically. The semi classical picture of a Neel-like ground 
state is completely off. This is special to 1 dimension and actually arises as a general 
consequence of Coleman's theorem it states that in (1 + 1) dimensions, infrared 
divergences associated with Goldstone bosons will always wash out the classical value of 
the order parameter rendering spontaneous symmetry breaking of continuous symmetries 
impossible. One is therefore forced to look elsewhere in order to describe the low energy 
physics of the Heisenberg model. 

1.3 Non-Linear a (NLa) Model 



The most consistent continuum model derivable from Eqn. (LI) is the Nonlinear a model. 
In addition to being a continuum model (valid only in the long wave length limit), it is 

— * 

also based on a large s approximation. One introduces two fields: 4> — corresponding to 
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the Neel order parameter, and 1 — the uniform magnetization. The spin variable, S is 
defined in terms of these fields via 

S{x) = {-iys${x)+\{x) (1.18) 

The conventional derivation defines these variables on the lattice and a continuum limit 
is taken in the semi-classical large s approximation to arrive at the NLcx model Hamil- 
tonian 0. There are several problems with this approach. First, parity is broken in 
intermediate steps and is eventually restored in the continuum limit. Second, and more 
importantly, the crucial topological term which is found in the continuum Hamiltonian 
is derived without clear notions of how 1/s corrections may be made. A much more 
elegant approach will be reviewed here. We will make use of path integral formalism 
based on independent derivations by Haldane [0 and Fradkin and Stone [^.^ These were 
motivated by similar questions about topological terms in 2-D quantum models thought 
useful in attempting to describe high temperature superconductors. 

One begins by defining a coherent basis of states for the S = s representation of 
SU{2) i: 

|n >= e*(^>^V 5 > (1-19) 

where | + s > is the eigenstate of with eigenvalue s, z ■ n = cos 6 and {z x n) is a unit 
vector perpendicular to both h and z. We see that \n > is the state with spin s in the 
n direction (ie. n ■ S\n >= s\n >.) This basis is over-complete and not orthogonal. We 
now make use of two identities: 

< n,\n, >= e-*(^-^^) ( il^ll^V (1.20) 



2s + 1 

4:71 



2 

J3 



d^n S{n ■ n - l)\n >< n\ (1-21) 



^We note that the derivation in reference g is somewhat in error. We correct their mistakes using a 
similar derivation due to Ian Affleck (unpubhshed). 
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$(ni,n2) is the area enclosed by the geodesic triangle on the sphere with vertices at 
ni, 77,2 and the north pole. There is an ambiguity of Air in this definition, but this makes 
little difference when exponentiated since e^'^** = 1 for s integer or half-integer. The first 
identity is most easily proved by using the ||, | > representation of SU{2) while the 
second follows from the first. We can now use these states to write the partition function 
or path integral of the system: 



Tre-^^= lim / 



'^-^ 2s + 1 

n ^d^Hn) 

Ll=l 



N 



n 



-AtH\ 



^(Tm-l) > 



J m=2 

With At — > 0, we expand the exponential to order 0(At) to get 



(1.22) 



< h{Tm)\e ^■^^|n(r^_i) >^< h{Tm)\h{Tm-l) > - < n{Tm)\H\n{Tm) > At 



= e 



In the limit N ^ oo, the path integral can be written, 



Tre-^^ oc J [Dn{T)] 



(1.24) 



E 

m=2 



< n{Tj\H\n(rm) > Ar - s fi$(n(r^), n(r^-i) + ln( ^ + ^^""^ ^ ^^""^^^ ) 



The last term can be written to second order in At as 



sAt rP 



rp o sAt ri^ 

/ dTh{T) ■ ^,h{T) ^ —- / dT{drh{T)) 

Jo 4 Jo 



(1.25) 



This vanishes in taking the limit At 0. The sum over the phases is just the area 
enclosed by the vector n{T) as it traces its periodic path on the surface of the sphere. 
Parametrizing n as 



n — (sin 9 cos 0, sin 9 sin 0, cos 9) 



;i.26) 
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we can write 



-is J dA == —is J (i0(l — 008 6*) = —is J dt {1 — cost 



;i.27) 



The Hamiltonian, Eqn. ( |1.1|) , is a sum over a chain of spins. We must therefore 
extend the path integral to all sites. This is done by indexing each of the coherent states 
with a position label, x, and making the substitution 



:i.28) 



Note that 



— * 

< n(r, x)| S{x)\n{T,x) >= sh{T,x) 



;i.29) 



It is useful to write n{T,x) in terms of a staggered and a uniform part which are slowly 
varying in the limit of large s: 



n{T,x) = {-iy${T,x) + 1{t,x)/s 



;i.30) 



To leading order in 1/s and derivatives of the slowly varying fields, this produces the 
constraints 



— * — ► — > 

4>{t, x) ■ 4>{t, x) = 1 <P{t, x) ■ 1{t, a;) = 



;i.3i) 



Setting Ax = 1 (the fact that the fields vary slowly over this interval is justified a 
posteriori), we find that the leading contribution of the Hamiltonian to the action in the 
continuum limit is 



J dx dr {^(d^4)f + aT / s' 



;i.32) 



We add up the phase terms by combining them in pairs: 



-isA 



—IS 



dx [A 0(r,x+l)+ r(r,x + l)/s +A -${t,x)+ r(r, ) (1.33) 
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Because A is an oriented area with respect to its argument, changing the sign of the 
argument also changes the sign of the area. Eqn. ( |1.33| ) can be written, 



-isA=^Jdx (a 4>{t,x) + 5${t,x) - A ${t,x) ) (1.34) 

— * — * — * 

where to leading order, S(I){t,x) = dx(f>{T,x) + 2 l/s. This then gives, 

— is f ~* / ^ ^ 
—isA = J dx d<p ■ (^0(r, x) x 50(r, x) 



—IS 



J dx dr 4>{t, x) ■ {6(f){T, x) x dr4>ij^ 



X 



J dx dr i 4){t, x) ■ {dx4>{T, x) x dr4>{T, x)) — 2— ■ (</>(r, x) x dr4>{T-, x)) J (1.35) 



2 



If we compactify cf) so that (f) constant for + r^| — >■ oo, and maintain the constraint 

-»2 

0=1 (valid to 1/s^), one can recognize the integral 

Q = 1-J dtdx$- (dr4) X dx4)) (1.36) 

as measuring the winding number of the sphere onto the sphere. The integrand is the 
Jacobian for the change of variables from compactified coordinates on the plane to those 

— * 

on 0-space (also a sphere). Q is an integer corresponding to one of the countably many 
topologically inequivalent ways there are to smoothly map the sphere onto the sphere; 
thus the phase term can be written as —2'KisQ. For s an integer, a sum over all possible 
topological configurations will not affect the path integral. For s half-integer, however, 
we can expect a drastic difference, since the path integral will be the difference between 
partition functions with even and odd Q's. It is important to stress that this is a purely 
quantum mechanical result which has no analogue in the 2-D finite temperature classical 
Heisenberg model (there is a well known equivalence between (rf, l)-dimensional quantum 
field theory and d + 1-dimensional finite temperature classical statistical mechanics |T(1|). 
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A detailed discussion of how a half-integer s will affect the physics will be omitted here; 
the reader is instead referred to and references therein. 

— * 

We can solve the equations of motion for I: 

1 = -1- (0 X dr4>) (1.37) 

Not surprisingly, I is the generator of rotations. After integrating out the I fields, the 
final action is, 

S = -2tiisQ + — dxdr [d^$) + ^ / {dvr$) (1-38) 
Where we now define, 

v = 2Js g = -, (1.39) 
s 

The action can be written 

S = 2msQ + — [ dxdr d^4>d^4) (1.40) 

It is clear how 1/s corrections entered into the calculation of the topological term. More- 
over, we did not break parity to derive ( [1.4(]| ). 

We are interested in integer s (in fact, s = 1). To this end we may ignore the 
topological term in the action, as discussed, and consider the nonlinear cr-model: 



£ = —d^4>d^'4> 4> = 1 (1.41) 



We now motivate the idea that, contrary to spin- wave theory, this model features a gap 
in its low energy spectrum. We first do this in the spirit of reference [0. We can deal 
with the constraint by incorporating it into the path integral as a Lagrange multiplier: 



Z oc / V(f)VXe V J (1.42) 
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The constraint is now enforced by the equation of motion for A. The (j) fields can be 
integrated out in the usual way to give 

ZocJ pAe^/'^'"^(*)-^i°s^^*ll-^'+^ll (1.43) 

where is the number of components of 0. As — > oo the path integral is concentrated 
near the smallest value of the argument of the exponential. Minimizing this argument 
with respect to A we solve for the saddle point, A: 

= <x\{ =- ) \x > (1.44) 

using the standard rules for functional differentiation. The RHS of the above is simply 



the Green's function for a boson field with mass vyX 

, f 1 1 , f d^k 1 , , 

< X\ < }\X >= / 7 TTT (1-45) 

'l-a2 + Aj' J {27ry kt^k^ + X ^ ' 

1 , A2 

= — log — 

47r ^ A 

where A is an ultraviolet cutoff, d^k = dk duo/v, and k^k^ = k"^ + u'^/v'^. Solving for the 
square of the mass. A: 

Another way to see the presence of a mass gap is to integrate out ultraviolet modes 
and apply the renormalization group. We start with the Lagrangian Eqn. ( |1.41| ) and 



parametrize the fluctuations in terms of slow and fast modes. One then integrates out the 
fast fields. This calculation is logarithmically infrared divergent. One then renormalizes 
by subtracting out the offending terms from the effective Lagrangian. Equivalently, one 
can achieve the same effect to the same order in perturbation theory by redefining the 
coupling constant in terms of its bare value. A calculation of this sort (for the 0{N) 
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model) is done in reference |jTO[- The renormalized coupling constant becomes 



ZTT 

With qq = 2/s, we now see that the coupling constant is of order unity for length scales 

e ^ e""' (1.48) 

Keeping in mind that this is a 'Lorentz invariant' theory, there must be a corresponding 
mass scale, A: 

A oc e-""' (1.49) 

There are other similar heuristic calculations that suggest a mass gap; none are iron- 
clad, but the sum of them together makes for strong evidence that indeed the s = 1 1-D 
Heisenberg antiferromagnet is disordered at all temperatures and is well described by 
the NLo" model. Better justification comes from exact S'-matrix results and numerical 
simulations. The exact S'-matrix results are due to work by Zamolodchikov and Zamolod- 
chikov |jTl|, and Karowski and Weisz |]12|. The 0(3) invariance of the NLa model allows 



for an infinite number of conservation laws. These imply strong constraints on S-matrix 
elements and, consequently, on on-shell Green's functions. One characteristic of such an 
S-matrix is factorizability. This means that A^-particle scattering can be expressed as 
products of 2-particle scattering matrix elements. The simplest such S'-matrix consistent 
with the symmetries of the NLo" model has a triplet of massive soliton states with an 
effective repulsive local interaction. This conjecture has been checked in perturbation 
theory in 1/A^ (for the 0{N) NLa model 0) to order l/N'^. 



Numerical results have been pursued since Haldane made his conjecture in 1983 |I3^ 



14| , |T5|, |T6[. They have all essentially confirmed Haldane's picture and the validity of 



the NLo" model. To date, the best numerical work has been due to White's method of 
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the density matrix renormalization group [|14| and recent exact diagonalization The 
former predicts a gap A = .41050(2) J, while the latter has A = .41049(2) J. Numerical 
investigations of the spin operator structure factor |T6|, S{k), for the isotropic chain 
show remarkable agreement with the 'exact' S-matrix result for two magnon production 
over a region larger than expected (two magnon production is known to dominate at low 
momenta, k < .Svr, from numerical studies |TB]. This can be probed in neutron scattering 



experiments fl^-) For higher momenta, one must include one magnon contributions 
which dominate as A; ^ tt. The intermediate region in momentum space, .Svr < k < .Stt, 

— * 

is not expected to be well represented by the NLcr model; this is because the fields, (p 

— * 

and I describe low energy (and therefore large wavelength) excitations about k = n and 
k = 0, respectively. The same study also determined the correlation length, ^ = 6.03(1), 
the velocity, f/J ~ 2.5 and the coupling constant, g ~ 1.28 in rough agreement with the 
1/s expansion result |l^ g ~ 1.44 and the value derived above, g = 2/ s = 2. 

This ends the introductory discussion of the NLo" model. A more in-depth approach 
will be taken when we consider anisotropics and develop the necessary tools to calculate 
the NMR relaxation rate in Chapter 2. 



1.4 Boson Model 

Although the NLa model is convincingly accurate in describing the low energy physics 
of the Heisenberg 1-D antiferromagnet, it has several deficiencies. First, off-shell Green's 
functions are not known; and second, anisotropics are not easily tractable within the 
framework of the model (the S- matrix is no longer factorizable, as earlier discussed, since 
one loses the infinite number of conservation laws). A happy compromise which contains 
all of the qualitative aspects of the NLo" model and yet allows for more computability 
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and generalization is the Landau-Ginsburg boson model pQ 



—2 

where the constraint 0=1 has been relaxed in the Lagrangian of the NLo" model, and a 
0^ interaction has been added for stability. The Hamiltonian, (|1.50|) , possesses the correct 
symmetries, three massive low energy excitations, and a repulsive weak interaction. As 
with the NLo" model, the field (p acts on the ground state to produce the triplet of massive 
excitations or magnons. We note that this model becomes exact in taking the N —>■ oo 
limit of the 0{N) NLcr model ||10| (recall that is the number of components of the 
field (j)). As in Eqn. ( |1.37|) , the generator of rotational symmetry (the uniform part of 



the spin operator) is 

l = 4>xn (1.51) 



where II = (we absorbed the coupling constant, (yf, into the definition of 4> in (|1.50|) . 



Expanding (f) in terms of creation and annihilation operators, we see that I acts as a two 
magnon operator producing or annihilating a pair, or else flipping the polarization of a 
single magnon. This picture is obvious in this simpler model, whereas the same analysis is 
only confirmed by the exact S'-matrix results and the gratifying agreement with numerical 
work in the case of the NLo" model. The gap, A can be phenomenologically fitted to 
experiments such as neutron scattering as can be the 'velocity of light', v. Including 
on-site anisotropy 

^aniso = E {D{Stf + E{{S^f - {8^')) (1-52) 

i 

simply amounts to introducing three phenomenological masses. This will be discussed in 
more detail in Chapter 2. 
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On comparison of the predictions of both models one finds overall qualitative agree- 
ment in studies of form factors |16|, 0. As one moves away from zero wave vector the 
agreement between the models weakens. This makes for one of the disadvantages of the 
bosons. Also, in attempting to calculate certain Green's functions, such as the staggered 
field correlation function, one is forced to rely on perturbation theory in A. Although 
A can be phenomenologically fitted, there is much ambiguity in choosing the interac- 
tion term. One can equally put in by hand any positive polynomial term in . This 
is because the fields carry no mass dimension making all polynomial interaction terms 
relevant. It should be understood that this model is phenomenological and is introduced 
for its simplicity. In the final analysis, justification for its use must come from numerical 
and real experiments. 

1.5 Fermion Model 

Before introducing the next model, we would like to begin by apologizing for the cryptic 
description of the concepts to be mentioned in this section. A deeper understanding 
would require a diversion into conformal field theory tangential to the main lines of the 
thesis. Instead, the reader is invited to investigate the literature. 

There is another model exhibiting some of the desirable properties of the boson model. 
This is an analogue of the Landau-Ginsburg model but phrased in terms of a triplet of 
relativistic fermions: 




1 - 



d 



1 -t 



d 



^^{^R >L - V'l • ^r) + K^L X ^l) ■ ii'R X ^r) 



(1.53) 
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The fields if) are Majorana (Hermitean) fermions with equal time anticommutation rela- 
tions 

= 5ss'S'^ Six-y) S,S' = L,R (1.54) 

The L and R label left and right moving fields, respectively. This model is not trivially 
related to either of the models described above; it was first introduced by Tsvelik |^ 
to achieve better agreement with experimental data on the anisotropic Haldane Gap 
material NENP. The motivation comes from a model sitting on the boundary between 



the Haldane phase and a spontaneously dimerized phase |22], with the Hamiltonian 



H = JY,[S,- - ( 5, ■ (1.55) 



This Bethe Ansatz integrable Hamiltonian features a gapless spectrum and has a con- 
tinuum limit equivalent to a A; = 2 Wess-Zumino-Witten (WZW) NLa model. This, in 
turn, is a conformal field theory equivalent to three decoupled critical Ising models. 
The well known mapping of the critical Ising model to a massless free Majorana fermion 
24^ brings us to write ( |1.55| ) as 



Reducing the biquadratic coupling in ( |1.55| ) moves the Ising models away from their 



critical point. Symmetry allows the addition of interactions corresponding to mass and 
four fermi terms, as in Eqn. (|1.53|) . The four fermi term proportional to A is the 
only marginal one allowed by 0(3) symmetry. It will generally be ignored or treated 
perturbatively, in a similar phenomenological spirit to that of the Landau-Ginsburg boson 
model (ultimate justification for this, as for the boson model, comes from numerical and 
real experiments). For weak interactions (which is the case assumed) all Green's functions 
will have simple poles at the phenomenological masses and will be trivial on-shell. The 
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off-shell behaviour depends on the interaction terms chosen and is therefore very much 
model dependent. 

It can easily be checked that (|1.54|) gives the right commutation relations for the 
SU{2) algebra, [l'{x),P{y)] = iS{x - y)e'^H''{x), with 

^=^(V5LX^L + V^i?, x^ii) (1-57) 



2 

— * 

This allows us to identify / with the generator of global rotations or the uniform part of 

— * — * — * 

the spin, S. Expanding ■i/?^ and i/'l terms of creation and annihilation operators, we 

— * 

see that, here too, / is quadratic in such operators. Notice that this representation for 

— * 

I does not couple left and right movers. This is in sharp contrast to the boson or NLcx 
models (where one can write the boson operator as a sum of left and right moving parts). 
We will later see that this point can potentially give experimental predictions which will 
contrast between the models. 

^ — * ^ 

The particles created by the fields, ip, are identified with massive magnons. The 
masses can be fixed by hand to agree with the experimental dispersions so that on- 



site anisotropy terms coming from Eqn. (|1.52|) can be easily parametrized, as in the 



boson model. Other interaction terms which might arise from breaking the symmetry 
are usualy ignored for ease of calculation. As always, ultimate justification for this is 
found in numerical and real experiment. 

As mentioned above, / is again a two magnon operator. It is also possible to represent 
the staggered magnetization (the analogue of (p) in this approach, but it is considerably 
more complicated (one can use bosonization techniques p3|). Near the massless point, 
this operator reduces to the fundamental field of the WZW model, or equivalently to 
products of the order and disorder fields, a and /x, of the three Ising models 
These operators are highly non-local with respect to the fermion fields. The corresponding 
correlation functions can be expressed in terms of products of Painleve functions 
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27| . They exhibit poles at the fermion masses together with additional structure at 
higher energy. Unlike the free boson model, a simple interpretation of the staggered 
magnetization density as a single magnon operator doesn't hold. This complicates the 
use of this model. 

One way to justify the use of the fermion model without resorting to complicated 
explanations is to notice that in the long wave length limit of the 0(3) symmetric case, all 
models are in agreement (see Chapter 2). For smaller wave lengths, the different models 
correspond to different continuum representations of the lattice model. 0(3) symmetry is 
broken differently in each model (for example, see the different results for matrix elements 
of I in Chapter 2). The idea is that we have three (two, for lesser symmetry) workable 
descriptions whose ultimate merits can only be decided phenomenologically. 



1.6 Nuclear Magnetic Relaxation Rate 

Experiments on condensed matter systems typically measure observables which are di- 
rectly related to Green's functions. This is no surprise since most such experiments 
measure the response of the system to an external probe. This is in contrast with par- 
ticle physics experiments which usualy examine the nature of scattering into asymptotic 
states. Formally the difference is that particle physicists measure time-ordered Green's 
functions while their friends in condensed matter physics measure retarded Green's func- 
tions. The nuclear magnetic relaxation rate, 1/Ti, measures the local correlations of the 
system at low frequency. The probe is the nucleus of some atom in the sample which 
has a non-zero nuclear magnetic moment weakly coupled to the system of interest. In 
the case of the Heisenberg 1-D antiferromagnet, we assume that in addition to the spin 

— * 

Hamiltonian, Hs, there is also Zeeman coupling, Hz, to a uniform magnetic field, H, 
by both the nuclear and Heisenberg spins, and that there is a hyperfine coupling between 
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the two systems, -f^Hyper- ^^^^ assume for simplicity that the nuclear spins do not 
directly couple to each other. 

^Total ^Hs + Hz + ^Hyper = + %yper 
= Hs-J2l^B H-Ge- S^-J2l^N H-G^- + Si-A,,- Ij (1.58) 

i 3 i,3 

Ge and Gat are the gyromagnetic tensors for the electron and nuclear spins, respectively. 
Aij is the hyperfine tensor coupling the nuclear spin on site j to the electronic spin on 
site i. We now define the characteristic frequency 

lun = A*jv| H\ (1.59) 

In nuclear magnetic resonance (NMR) experiments one strives to temporarily achieve 
a non-equilibrium population difference between nuclear spins with different spin eigen- 
values along the uniform field direction. This is normally achieved with pulses of RF 
electromagnetic radiation possessing ac magnetic fields perpendicular to the externally 
applied uniform magnetic field. As is well known, a resonance phenomenon occurs at 
RF-frequencies near ujn (in reality it is easier to tune the uniform field to resonate with 
a fixed RF field). 

In the presence of a non-equilibrium occupation of states, the nuclear spins "relax" 
towards an equilibrium configuration by making transitions between states of different 
spin eigenvalues. This would not normally be possible if the nuclear spins were completely 
free. Coupling to another system is necessary in order to conserve energy during the 
transitions. The energy given off or absorbed must induce a corresponding transition into 
a different energy state in the system which couples to the nuclear spins. Let us illustrate 
the situation with an s = 1/2 nuclear spin. In the absence of hyperfine interactions, we 
assume that is a good quantum number (where z is the direction of the static magnetic 
field), and that Gn is isotropic (these are generally good assumptions). The rate equation 
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for the number of nucleii, A^_|_, with = +| is 

= -N+n+^_ + N_n_^+ (1.60) 

Where the transition probabihties per unit time are given by fij-^zp. We can rewrite this 
in terms of the total number of spins, N and the population difference, n: 

d Ti 

— = N{Q^^+ - n+^-) - ^^(^]_^+ + (1.61) 

Now, if we define 

^ ^ (fi_+ + Q+^_) (1.62) 
1 

then we see that in the limit that the transition rates only depend on time scales much 
shorter than those characteristic of the experimental probe, and in the limit of linear 
response (ie. Q is independent of n) the solution to (|1.61|) is 

n(t) = no + a e-*/^i (1.63) 

Where Uq is the equilibrium population difference (at finite temperature, states with an 
energy difference will necessarily have a population difference). We see that the relaxation 
rate, 1/Ti, describes the evolution of the nuclear system towards thermal equilibrium, or 
likewise, the decay of the population inversion magnetization achieved by RF pulses in 
NMR. 

We now derive an expression for the rate, 1/Ti. For a system with more general 
/, 1/Ti for a transition from an initial state with P = m to one with = m + 1 is 
normalized by the factor /(/ + 1) — m{m + 1). To begin, we need an expression for 
the transition rate Qm->m+i describing a nuclear spin at site j starting in the state with 
P = m and ending up with P = m + 1. It does not matter which nuclear spin we pick 
if we assume translational invariance; since there is no nuclear spin-spin coupling'^, the 



■^In NENP, the dipolar nuclear spin-spin couplings are roughly 200 times smaller than the hyperfine 
coupling. 
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relaxation rate for one is the relaxation rate for the whole system. Let us assume that 
the initial and final Heisenberg spin state are given by the labels n and n', respectively. 
Then Fermi's Golden Rule gives 

^7?=m,n— »/?=(m+l),n' 

27r| < /| = (m + 1), n'|i7fjypgj.|i"| ^m,n> - — (1.64) 

Notice that we multiplied the normal expression for the Golden Rule by the Boltzmann 
probability that the Heisenberg spin system is in the initial state, \n >. The only part 
of -f^Hyper '^^ich will contribute is S^A!^j~Ij': 

^/|=m,n^/| = (m+l),n' — 

27r(7(7 + 1) - m(m + 1))| ^ A^,-- < v!\S^,\n > - K - c^iv)^^ (1.65) 

The analogous expression for il(m+i)n'-»/?=m,n is 

^7?=(m+l),n'^7?=m,n — 

2tt{I{I + 1) - m(m + 1))| ^ A^/ < n\S:^\n' > \^S{En' - - un)^^- (1-66) 

1,1/ 

Since A^^"*" — (A^^")* we get for the relaxation rate 

(Vri)„w = 'lAHKf < n'\Snn > mE^^ - K - c^iv) ^' + ^ " ) 

We now sum over all possible transitions to arrive at 1/Ti: 

^ = 27rE I E Ar < ^'l-^ri^ > \'S{En' - E^-u^)^- ^ ^ (1.68) 

1 n,n' i,!^ 
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Since the sum over the states n and n' is a trace over states in the Heisenberg spin system, 
we can conveniently restrict ourselves to that system only and write 

W = r '^^ < [ll^'tfS-itlY.^'rnPm] > (1-69) 



where <> denotes a thermal average. This is the famous expression derived in As 



promised, when Aij is well localized, we see that 1/Ti is related to the low frequency 
local correlation function. 



Chapter 2 
Details of the Models 



In this chapter we discuss in detail the three models introduced in the last chapter. We 
will derive the necessary tools to calculate the relaxation rate 1/Ti and mention some 
pertinent issues which can be important in investigating 1-D Heisenberg antiferromagnets 
(IDHAF's) using other means. 

2.1 NLcT Model: Temperature and Field Dependence of the Spectrum; Exact 
Results 

2.1.1 Temperature Dependence of the Gap 

In this section and the next we will discuss how the excitation energies of the lowest 
modes change with varying parameters. This is especially important when one chooses 
to perform calculations using the Landau-Ginzburg boson model. Since this model im- 
plicitly adopts the gap parameters from the NLcr model, any dependence of the gaps on 
magnetic field or temperature must first be calculated within the framework of the NLcx 
model. The results can become useful in interpreting experimental data using the simple 
boson model. 



We would like to begin by extending some recent work by Jolicoeur and Golinelli |^| 
on the temperature dependence of the low energy spectrum. It may seem strange or 
even contradictory at first sight to speak of a spectrum as being temperature dependent. 
What one must keep in mind is that the low energy description of the NLo" model as three 
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massive bosons with relativistic dispersion, is an effective one. The true excitations of 
the model are collective, and if we insist on maintaining a single particle description we 
should not be surprised that the effective single particle interactions will be temperature 
dependent (as, therefore, will be the effective one-particle spectrum). A similar approach 
is taken in BCS theory where the BCS gap has a temperature dependence arising from 
a consistency condition. 



In Chapter One we introduced a consistency equation, Eqn. ( |1.44| ), for the classical 
or saddle point value of the Lagrange multiplier field. A, in the NLcr model. This result 
always holds to lowest order in the fluctuating field A(x), regardless of the value of N in 
the large N expansion. Of course, it only becomes exact for — > oo. We can also look 

— * 

at the consistency equation as a constraint equation guaranteeing that the <p two point 
function is unity when evaluated at the origin; when = 3, 



l=<^W.^W>=G» = 5;/||j^ (2.1 



where we've assumed a renormalized mean value for A fT^. Notice, also, that we're 
choosing to work in Euclidean space. One can likewise see that the constraint equation 
is nothing more than a minimization of the zero point energy of the system with respect 
to the fluctuating field A: 

dk 



Eo = 3 J—Uk - A J (2.2) 



where Uk = v^k^^k^ + A. If we choose to add on-site anisotropics to the model, as 
in Eqn. ( |1.52| ), then the contribution to the Lagrangian (modulo irrelevant terms which 
also break 'Lorentz invariance') is 



-D{Stf - E{{S:f - {Sff) - -D{<p^{x)r - E{{r{x)f - {<P\x)f)) (2.3) 
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We can always find some axes, xyz, so that the addition to the Lagrangian is in the 
above form. We can read off the new constraint equation on the Green's function 

^ v"^ f dk dtu f 1 



J J (27r)2 I + v^k^ + A2 + Al 



D 



Where we've once more assumed renormahzed masses with the correspondence, v'^X + 
2vD ^ A^ + A|, and v^X ± 2vE A^ ± A|;. As usual, making this model temper- 
ature dependent consists of replacing the integral over u; by a corresponding sum over 
Matsubara frequencies. Eqn. (|2.4|) becomes. 



_Tv^ r dk f 1 

~ 7 2^ / 



J 2n \ujl + v^k^ + A\T) + Al 



w2+i;2p + A2(T) + A| Cu2+t;2A;2 + A2(T) - A|^ 

In summing over the frequencies we use, 

1 ^ 1 l cot(^; 

^ ^ u;2 + f 2A;2 + m2 ~ ~2 iu;fc 



' + (2.6) 



2uJk \ el^'^k _ 1 



where uj^ = yv^K^^^^rri? . We also need the following two integrals 

1 / 2 \ r p-fSm-f3k^/2m 

dk — ]^2 dk , 



f V Tn 



r dk 9 

— ^- \og{2Kv/m) {21 

J UJh V 
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where A is an ultraviolet cutoff and we made the approximation /5m ^ f . Gathering all 
of the above, we can set Eqn. ( p.4|) equal to Eqn. ( p75|) to arrive at 

-A-/T1 

(2.9) 



, /A2(T)Al(T)A2(r)\ , " 



g-A,/T e-A+/T e-A-/T 



+ 



/a: 



+ 



/a: 



where A2 = y A^ + A|) and A± = y A^ ± A|;. Eqn. ( |2.9| ) implies a cubic equation for 
the square of the temperature dependent gap, A^(T). Once more making the approxi- 
mation, /3A_ ^ 1 we can linearize the equation and solve for A(T): 

-A+/T 



A(T) - A 



V27fr 



-A_/T 



/A^ 



+ 



+ 



In the 0(3) symmetric case this reduces to the formula derived in ||29|. We would like to 
say a few things about Eqn. ( p^.lUD before going on to the next section. First, notice that 
we implicitly assumed that only the expectation value of the fluctuating field A acquired 
a temperature dependence. The renormalized values of the anisotropies, Aj^ and A|; do 
not. This is because renormalization occurs at T = first. At non-zero temperatures the 
free energy may acquire a term linear in the fluctuations of A; the constraint equation, 
Eq. (|2.1| ), amounts to cancelling that contribution in the Lagrangian. The exponential 
terms logically appear as a result of calculating 



< A >= Tr 



(2.11) 



^2 



and then subtracting < \ > cp from the Lagrangian. We would also like to point 
out that the validity condition for this analysis, f3Am ^ 1, where Am is the smallest 
gap, is more robust than seems. It is well known that for a value of the anisotropy. 



-D ~ J 0, |T8[, the lower gap closes and the system goes through a critical point, into 
a phase with a new singlet ground state (the order parameter is a non-local operator 
in spins, and in fact, this transition is not reproduced correctly via the NLcr model) 
and a gap. At large negative values, D ~ —.4 J, the system goes through an Ising 



Chapter 2. Details of the Models 



29 



transition into an antiferromagnetically ordered phase. The bottom hne is that \D\ must 
be small in comparison with A. Similar, but more obvious, cautionary remarks apply 
to E. Moreover, as discussed in ||2^, the NLcr model is not expected to remain valid 



at temperatures of order twice the gap. This is because the model does not exhibit a 
maximum in the heat capacity and in the magnetic susceptibility as shown in numerical 
studies at these temperatures. 

Finally, notice that the difference in gaps will close as T increases. This is no surprise 
since at high temperatures the mass scales are irrelevant and we expect a restoration of 
0(3) symmetry. 

2.1.2 Field Dependence of the Gap 

We will be interested in adding a magnetic field term to the Hamiltonian. To do so 
consistently, we must couple the magnetic field to the generator of rotations (the total 

— * — * 

spin operator), J2i Si. In terms of the continuum fields, we couple the magnetic field, H, 

— * — * — * _ 

to I via, —Qe^-B H ■ J dx l{x), and add this to the NLcr Hamiltonian. The corresponding 
Euclidean space Lagrangian is 

C = -— (\d$/dt + /I X 0|2 + v\d4>/dxf - 2vD{(t)'f - 2vE{{(j)''Y - {^y)^)) (2.12) 



— * — * 

where h = Qe^-B H. In the 0(3) and U{1) symmetric cases, where the field really couples 



to a conserved charge, no other terms are allowed in ( |2.12| ). In the case of lower symmetry, 
we retain this as the simplest form, realizing that other, more complicated terms may 

— * — * 

arise. This time, when we integrate out the fields, the eigenvalues of the propagator 
are not as trivial. However, if we assume that the field is placed along a direction of 
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symmetry, say the 2;-direction, then the Tr log will be over eigenvalues of the matrix 
/ uj^ + v^e + A^-h^ + Al 2ujh ^ 



-2iuh cu^ + v^k^-h^ + A^-Al (2.13) 

y U U UJ~ + V"K." + ZA" + ZAJ, j 

Where, again, we've assumed renormalized values for the masses. The eigenvalues, r/j, 
are 

r]:i = u^ + v^e + A\h) + Al 



= u;^ + v^k' -h' + A\h) ± J-Ah^u^ + A% (2.14) 



We can write the constraint equation as 
cPk ( 1 



O. 0-»0. AO. AO 



2 



(27r)2 \^Cj2 + t;2A;2 + A2 + + t;2p + A2 + A| Cj2 + t;2A;2 + A2 - A^ 

d^k l\ \ \\ 
TTTT^ - + — + — (2.15) 

The integrand on the right hand side has poles at the negative solutions of the equations 
of motion 

u;l={e + A\h) + Al) 



uj^^=ie + h' + A\h) ± J4/i2(A;2 + A2(/?.)) + A% (2.16) 



Integrating over these poles gives 

fdk(j_ uol-k^ - A\h) + h^ - k^ - A\h) + h'' \ 

Before continuing, we mention that the T dependence can easily be worked in by mul- 
tiplying the terms in the integrand above by 1 + -g^S^:-^, respectively. It is possible to 
simplify Eqn. ( |2.17D further to read 



dk ( I I I Ah^ \ , , 

^ - + — + . ^ J (2.18) 
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where it is now clear that ( p.l5| ) is satisfied for h ^ 0. This can be shown to be the same 



result obtained by minimizing the zero point energy. To solve for A{h), one must decide 
on a sufficiently large ultraviolet cutoff and resort to a numerical root finding routine. 
In Figure |2.1| we plot the energy gaps of Eqn. ( p. 161) unconstrained by Eqn. (|2.15|) . To 



compare, we also plot the field dependent gaps of the fermion model which are generally 
considered in agreement with experiment |^ (at least for the material NENP). The zero 
field gaps are fitted to the gaps found in NENP: = 1.17meV, A;, = 2.52meV and 
Ac = 1.34meV. The subscripts, 'a, 6' and 'c' refer to appropriate crystal axes of NENP. 
The dispersions differ most at higher fields, and for large A^;. In Figure p.l.2| we replot 



the gaps but this time correct for the constraint implied by Eqn. ( |2.17| ). The lower 



branches show good agreement right up to fields close to critical. There is, however, a 
greater discrepancy in the gap corresponding to the field direction. 

We now turn our attention to a seeming infrared catastrophe which occurs as h 
approaches the critical value given by 



/i, = ^A2(/i,)-A|; (2.19) 

this is where the lower gap closes and the integral ( |2.18| ) diverges logarithmically in the 
thermodynamic limit. At first sight one may hope that for k = 0, the last two terms 
in (|2.18|) conspire to eliminate the divergence for some value of h and A(/i) satisfying 
( p.l9| ). This, however, requires that 



5 4 



9Ai 



(2.20) 



9 9^ 16A^{K 

be simultaneously satisfied; this is impossible unless E = 0. In fact, for E = 0, A{h) 
is independent of h. This can be seen directly from Eqn. (|2.18| ) or by understanding 



that the variation of the zero point energies uj± = y/v'^k'^ + A^ ± h with respect to A^ 
is independent of h. Let's try to get a deeper feeling as to what's happening. Instead 
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Figure 2.1: Free boson and free fermion dispersions with the gap parameters of NENP. 
Top graph: H\\b; bottom graph: H\\c. 
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of starting with a finite E we can place the magnetic field near the U{1) critical value, 
he — A(/ic), and turn on the anisotropy. We write 



We now see that the limits h ^ he and E ^ don't commute after taking the derivative, 
Regardless of the limit at which we start, we can expect trouble when 



This is true even for /i = 0, which is the simplest case where this problem appears. 

Essentially the trouble arises because, returning to the NLcr action, we integrated out 
low energy modes. In the U{1) case the day was saved by symmetry which prevented the 
variation of < A > with magnetic field. There is no such luck in the Z2 x Z2 x Z2 scenario 
and for a sufficiently large value of it is no longer correct to integrate out all the fields 
even for large N . In the case of large N , one can integrate out all modes but the gapless 
one to arrive at a critical theory. This should be done when log (^-^^^^ ~ N/g. For 
N — 3 this criteria may be too restrictive, and instead one can adhere to £■ ^ A(/i) — h 
as the regime where the analysis of this section is valid. 

2.1.3 Exact Results 

In this section we note some important results which will be useful in the calculations 
of Chapter 3. As was mentioned in the last chapter, the 0(3) model possesses special 
properties that allow for some integrability. In particular, at long wavelengths one can 
say much about the matrix elements of the spin operator even after the 0(3) symmetry 
is broken to C/(l) by a magnetic field. 

As will become clear in the next chapter, wc are mostly concerned with matrix el- 
ements < k,a\S\0)\q,b >, where < k,a\ denotes a single magnon state created by the 



u;^_{k^O) = \{A{h)-hf 




(2.21) 



4hA 



Er.A_{h)-h 



(2.22) 
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staggered magnetization operator, (pl (this acts as a free boson operator in the large N 
hmit) with norm 

< k,a\q,b>=2TT5ab5ik- q) (2.23) 

a and b denote polarizations of the magnon states. Clearly, < k + q,a\(f)^{0)\k,b >= 0. 
This is obvious when the bosons are completely free; interactions do not change this 
picture since they must all be even in bosonic operators. The two magnon operator, 
/*(0), is expected to contribute.^ The matrix element is given by the Karowski and Weisz 
ansatz [0, 

< k,a\r{0)\q,b>= le'^' ^'^'^^'^'^ Gie) (2.24) 



where Uk = \/k'^v'^ + and the rapidity variable, 0, is defined via 

sinh(^' — ivr) = k 
sinh(^" — iir) = q 

e = e' - e" (2.25) 

with 



dx (e-2^' - 1) sin^[a;(i7r - 6)/ 2 



This ansatz is believed to be exact for the 0(3) NLa model, but is only approximately 
true for the s = 1 Heisenberg model; however, numerical simulations |]16| are in excellent 



agreement with this form at least half way through the Brillouin zone. Since we will 
largely be interested in \k — q\v < A ^ ttv, this 'exact' expression is more than sufficient. 



^Strictly speaking, since I is quadratic in 0, we expect it to be a two magnon operator only in the 
large N limit. 
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Some particular limits of interest are k q and k ~ —q, corresponding to forward 
and back scattering, respectively; in the former case, 6 ^ m while in the latter, 6 ~ 
i7r + 2A;^;/A. 

G(i7r + 2 W A) ^ 1 - + ^) (5) ' (2-27) 

This expression is a different result than for free bosons and reflects the effects of inter- 
actions. We will later see how this might affect experimental predictions. 

The above results change when a D type anisotropy is added. Essentially, the dif- 
ference is that the function, G{0) changes and the gap in the energy of states on the 
5"^ = branch will be shifted. There are no exact results for this model which is why the 
phenomenological boson and fermion models are important. We can, however, give some 
universal (ie. model independent) results which only depend on the conservation of total 
spin in the z-direction, when the momentum exchange is small, v\q — k\ <^ A. The only 
part of 5'^(0) which will contribute will be, essentially, S^_^ ^ J dx S^{x). Since this is 
a conserved operator, we can write down the one-particle matrix elements immediately: 

< k,s^ = ±l\S^{0)\q,s^ = ±l >Ri ±1 (2.28) 

Also true for all values oi v\q — k\ are the following 

< k,s'^0\S'{0)\q,s'^a >= 
< k,s'^+l\S'{0)\q,s' = -l>=0 (2.29) 

In Chapter 3 we will show that it is largely this universal behaviour which determines 
the relaxation rate, 1/71 in anisotropic media. 

Note that adding a magnetic field to the 0(3) or U (1) system will break the symmetry 
in the 0(3) case (we naturally assume that in the U{1) case, the field is placed parallel to 
the U{1) axis), but will hardly change any other results in either model. This is because 
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in adding a magnetic field all we have done is add a term to the Hamiltonian proportional 
to the conserved charge Zg = J dx l^{x). Since Zg commutes with the Hamiltonian, they 
can be simultaneously diagonalized with all low lying states labeled by their Zg quantum 
number, +1,0 or —1. Matrix elements of operators can only differ, at most, by some 
cumulative phase which corresponds to turning on the field sometime in the past. 

2.2 The Free Boson and Fermion Models 

We now turn our attention to the phenomenological models introduced in Chapter 1. 
After introducing the formalism which we will require to perform calculations, we will 
compare and contrast the energy spectra and fundamental matrix elements. 

2.2.1 Diagonalization 

To do the necessary calculations we need to have a basis of eigenstates for the non- 
interacting Hamiltonian and know the expansion of the field operators in terms of cre- 
ation/annihilation operators for these states. 



For now we assume that the mass and gyromagnetic tensors, D and G respectively, 
are simultaneously diagonalizable and work in this diagonal basis (this is rigorously true 
when the crystal field symmetry is no lower than orthorhombic - a sketch of a proof is 



Bosons 



-.2 

Relaxing the constraint 0=1 in (|2.12|) we see that we seek to diagonalize 




(2.30) 
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found on p. 750 of gj). 





( 


A? 








\ 




/ 


9i 





\ 


D = 














G = 







92 







\ 








A^ 


) 




v 








^3 / 



Also, we mention that we've set = 1 = a, where a is the lattice spacing. This has the 
effect of measuring energy in units of inverse seconds or inverse mass: 

[E] ~ [s]-i ~ [M]-^ ~ [v] ~ [0]^^ (2.32) 

Diagonalizing ( p^.30| ) is tedious (especially when the field does not lie in a direction of 
symmetry, for then all the branches mix) but the idea is to find the right Bogoliubov 
transformation. Working in momentum space, we define 



4>{k,t = 0) = -j=[ aU + Sk] (2.33) 

n{k,t = 0)=t^[aU- a,] (2.34) 

[ai,ai]]=27T5ijS{k-k') (2.35) 

ujq is an arbitrary quantity with the dimensions of energy. We need such a quantity to 
represent the fields (f) and J7 in terms of creation and annihilation operators. It turns 

— * — » 

out that when one writes 0^ and Ilk in terms of the creation/annihilation operators 
which diagonalize the Hamiltonian, the dependence on ujq disappears. Furthermore, the 
eigenvalues of the Hamiltonian are also independent of ujq, as might be expected. We 
will now restrict ourselves to the case where the field lies in a direction of symmetry. 
This leaves ( |2.3CID with Z2 x Z2 symmetry. Now only the excitations transverse to the 
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direction of the applied field mix and we need only solve a (4 x 4) set of equations for 
the diagonalizing creation/annihilation operators. Without loss of generality, we take the 
field to lie in the 3-direction, and set = 1. All told, the Hamiltonian, Eqn. (|2.3CI| ), for 
the mixed states is 

/oo 
Hk dk 
-oo 

Hk = d\Adk + a^_kA*d_k + d\^d^_k + afcB*a_fc (2.36) 



A=^l+^-]-ha, (2.37) 
4 4u;o 2 



B = -^I + A (2.38) 
4 4u;o 

where cr2 is the usual Pauli matrix, and 

A? + v^k^ 
K = I I (2.39) 

/\l + v^k^) 

The momentum space Hamiltonian can be written in terms of a single matrix M: 

Hk = ( 4, a^k) Ml J (2.40) 



As discussed in we seek the eigenvectors of the (non-hermitian) matrix rjM, where 



A B \ / 1 

VM=\ \ "^=1 I (2-41) 

-B* -A* J y -1 

This comes from requiring the new diagonal creation/annihilation operators to have the 
standard commutation relations. This also imposes the unusual normalization condition 
on the eigenvectors: C V C = ^- 
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Summarizing the above, we need to solve: 



/ <£0T _|_ 

4 iujo 



-I 



\ 



t^o J _ K 

4 4aJo 



_l_ K 

4 4(^0 



-I 

2 



which can be manipulated to give 







{uJq — u)l — h<J2 —{yjQ + Uj)l — h(T2 
Ujl — hcT2 — + Uj1 + h(T2 



K 



(2.42) 



(2.43) 



The eigenvalues of riM are already known, as they are the solutions to the classical 
equations of motion and come in pairs ±uj±. These are naturally the same frequencies 
given in Eqn.(^]T^) with the proper substitutions made for the gaps 



U;l = k' + }l' + 



2 , A? + Ai 



± 



Furthermore, we need only work to find one eigenvector of each pair because if 



(2.44) 

^ X 
, Y 



is a right-eigenvector of r]M with eigenvalue cu, then | | is a right-eigenvector of 

X* 

rjM with eigenvalue —u (X and Y are themselves two-component vectors). 
The bottom rows of (|2.43|) give the following set of equations 



(2.45) 



= (^2 + v'^k^)x±,2 - (^o(^±^±,2 - ihuJo^±,i 



(2.46) 
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where it turns out to be convenient to work with x = X + F and ^ = X — Y . The top 



rows can be manipulated to give 

= (/i^ - a;|)x±,2 - ihu;o^±,i + u;oUJ±^±,2 (2.47) 

= -{h^ - ujI)x±,i - ihuJo^±^2 - ^o<^±C±,i (2.48) 
These can then be worked to give 

(A^ + v^k^ + u;l- h^)x±,i = 2u;o^±e±,i (2-49) 

(A2' + vV +u;l- h^)x±,2 = 2u;o^±e±,2 (2.50) 

(/i^ + + v'^k'' - ujI)x±,i = -2ihuJoi±,2 (2.51) 

{h^ + A^ + v'^k'^ - ujI)x±,2 = 2i/iu;oC±,i (2-52) 



Note that if we fix the phase of xi to be real then ^1 must also be real as X2 and ^2 must 
be pure imaginary. The normalization condition, X — Y^Y — 1, now allows us to solve 
for the eigenvectors which form the columns of the transformation matrix between the 
old and diagonal bases of creation/annihilation operators. In terms of x and ^ this is 

X±,i^±,i - X±,2C±,2 = 1 (2.53) 

The solution is 



I ^o^±(^' + ^1 + ^'k' - ^l) .2 54) 

' (A? + T;2A;2+a;2 _;,2)(;,2 + ^i+M + ^;2^2_^2)y 
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X±,2 



2ihujQ 



(2.56) 



^±,2 — -Z X±,2 



(2.57) 



With the inclusion of the trivially diagonal unmixed component (ie. the three com- 
ponent), we can define the three by three matrices x ^i^nd ^ with the columns labeled 
by the eigenvalues (+, — ,3) and the rows labeled by the original masses (1,2,3): Xsa = 
and = 



_aJO_ 



1 

au = — 

fc 2 



. One easily verifies that 



rt 



(2.58) 



Where the 6's are the operators which diagonalize H. Equations ( |2.54| )-( |2.58| ) are the 
main results of this section. Before ending, we give some limiting forms for x and ^. In 
the limit h 0, 



( 



x = C 



-If 



i 







v 
















(2.59) 



While in the hmit A2 Ai 

/ 

X = = 



V 



2VA2+«2fc2 
2VA2+i;2fc2 





2VA2+D2fc2 
2VA2+t,2fc2 





^A2+t,2A;2 y 






(2.60) 



^A2+l,2A;2 y 



Fermions 



We would now like to repeat the diagonalization procedure for the fermion model. The 

— * 

free Hamiltonian with minimal coupling to the magnetic field (ie. coupling only to I) 
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and simplest parametrization of the mass terms (corresponding to anisotropies and giving 



zero field dispersion branches cu^ = \/A^ + k^v^) is 



1 r 



i=l 



with, setting v = 1, 



'4'R 



°° dk 
2^ 



(2.61) 



(2.62) 



Jo Zrr ^ 



(2.63) 



{4, ait} = 2n6ij6{k - k') 



(2.64) 



Notice that we coupled the magnetic field to the generator of global rotations, 

— * _ 

/ dx l{x), given by ( |1.57| ). the Hamiltonian density in /c-space becomes 



Hk = dlMk dk 



(2.65) 



where 



M 



Ik — ihx iA 
—iA —Ik — ihx 



(2.66) 



(2.67) 



a 



L,fc 



The idea now is to diagonalize this matrix and find the eigenvalues and eigenvectors. 
In other words, find the unitary transformation which diagonalizes H. Once more, we 
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assume the field is in a direction of symmetry so tliat we need only diagonalize a 4 x 4 
matrix. Given that the field is in the 3 direction, the eigenvalues are: 



cu^ = k + h -\ 



± 



(2.68) 



It may be more illuminating to write out M in a basis that is more natural to the 
C/(l) problem. Using 



^R,k 




(2.69) 



2 t 



In this basis, M becomes 



M 



1 

V2 



1 1 



-I I 



( .+ t \ 



(2.70) 



(2.71) 



Where A 



Ai+Ao 



and b 



k\ — hcr3 iAcTi + iSl. 
^ — iAcTi — iSl. —kl + has 

. The equations for the components of the eigenvectors 



possess the symmetries 



1*1 <-> U2, Us <-> 1*4, /i <-> —h 



(2.72) 



Ui Us, U2 U4, LU —LU 



(2.73) 
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where uu is the eigenvalue. After some algebra, 

2/ A(/,- - ^0 , , 

2iS{k -uj) , , 

2zA(A; + a;) 

Using the normalization condition, 

4 

El^^r^l (2-77) 
1=1 

we set the phase of ui to be real for positive eigenvalues; the above symmetries allow us the 
freedom to choose a convenient phase for the Uis corresponding to negative eigenvalues. 

ui = 2d{k + uj){uj'' + A^ - {k + - S^) ^ (2.78) 
[AS'^ik + cuYiiuj'' + - (A; + - S'^f + AA''{k - cuf) + 

{{CU + hf -k'' + S^- A2)2((^2 + ^2 _ ^ _ ^2)2 ^ 4^2^^ ^ ^)2)|1 

We define the 6x6 diagonalizing matrix with columns given by the eigenvectors 

as 

= y^.^,^..,<-J (2-79) 



4 = C/X^ (3 (2.80) 
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U 



1 
72 



1 1 





i -i 

^2 

1 1 

-i i 

j 

The diagonal operators, are defined as: 



V 



/3a 



4 




(2.81) 



(2.82) 



Our freedom in choosing the phases for the eigenvectors corresponding to negative 
eigenvalues allows us to write 



X 



^ R T 



T R 



(2.83) 



Each index of this matrix runs over six states; the first and last three correspond 
to right and left movers respectively. In the case of U{1) symmetry or higher, each set 
would correspond to states of definite spin. 

The d's and c's correspond to left and right moving fermions, respectively. This 
becomes clear in the limit Ai — > A2 — > 0. Some limiting forms of R and T are: 
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T{h^O) 



1 uj-i — 


k 




1 LO\ — 


k 



LU2 



W2 








t<J3 / 



(2.85) 



R{5 0) 



1 

71 



LUA 













/ 



(2.86) 



T{S ^0) = — 



1 

V2 



\ 








U3A 












(2.87) 



2.2.2 Discussion: Comparison of Spectra and Spin Operator Matrix Ele- 
ments 



The spectra for the boson and fermion models are given by Eqns. (12^4^ ) and (12^681 ), 
respectively. In the case of U{1) symmetry, Ai = A2, the two sets of formulae agree. 
However, with the lower orthorhombic symmetry, the two models are in agreement only 
for low magnetic field, h <^ Min(Ai, A2). The difference is most significant at the critical 
field where the lower gap vanishes. The boson model predicts he = Min(Ai, A2), while 
the fermion model gives he = A/A1A2 (see Fig. and p.l.2|) . Experimental evidence 
seems to favour the fermion model, but there are some subtleties which have previously 
been ignored. The data supporting the fermion dispersion comes from neutron scattering 
and NMR relaxation rate experiments performed on the anisotropic IDHAF material, 
NENP. In analyzing the data, however, crucial structural properties were neglected in 
the interpretation (namely, the fact that the local chain axes did not coincide with the 
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crystallographic axes). This, we believe, also led to a seeming contradiction with ESR 
data on the same substance which seemed to side with the boson dispersion [0.^ Aside 
from material properties, the possible temperature and field dependence of the mass 
parameters, Aj, has also been ignored so far. Since the boson model derives from the 
NLcr model, one should incorporate such field and temperature dependence into these 
basic parameters. We showed that considering field dependent masses brought closer 
agreement on the lower branch dispersion between the models up to fields given by Eqn. 

We mention in passing that the Hamiltonian, Eqn. ( p.30| ), is not the only quadratic 
one possible when the magnetization density is no longer conserved. It is possible to 
construct a modified boson Hamiltonian including extra terms designed to reproduce 
gaps identical with the fermion model [Q. The only constraint on such terms is that 
they do not mix the = modes corresponding to the degree of freedom parallel to the 
field. It is not obvious, however, what justifies such a modification other than a more 
convenient spectrum which replicates the fermion model at low energies. 

The fermion model is expected to become more accurate close to the critical field. 
The nature of the critical point was established in Ref. |Q. With U{1) symmetry, the 
phase transition is in the two dimensional xy universality class. The lowest lying mode 
of the Landau- Ginsburg boson model can be reduced to a single free boson (a phase field 
corresponding to the Goldstone mode), but the parameters of the resulting low energy 
Lagrangian must be renormalized to give the correct critical exponents of the xy-model. 
One does not have to resort to such lengths with the fermion model which correctly 
describes the transition without interactions. This is expected on several grounds. First, 
the many body ground state wave-function for a dilute gas of repulsive bosons is simply 
that of free fermions multiplied by a sign function to correct for the statistics. Second, 



^For more details on these matters, please see sections 5.1 and 6.3. 



Chapter 2. Details of the Models 



49 



the U{1) fermionic modes can be represented as particles and holes using a single Dirac 
fermion with chemical potential h (this can be seen in the matrix equation (p.71| )); 
this means that at h > A the ground state will be occupied by fermion states, each 
with = 1, and hence non-zero magnetization. The simplicity of the coupling to h 
guarantees that interactions will be as important near he as they are near h = 0. In 
particular, they will be negligible in the dilute gas limit. We thus see that interactions 
become progressively more important close to criticality in a boson theory, while the 
opposite takes place in an equivalent fermion theory. 

In the Z2 X Z2 case we expect an Ising-like transition corresponding to the breaking of 
one of the Z2 symmetries remaining. Here things are even clearer. Mean field theory for 
the boson model is completely hopeless as is evidenced by the unphysical behaviour of the 
lowest lying gap at h > Min(Ai, A2). This function always possesses a zero even at non- 



vanishing fc-vectors. Moreover, it is imaginary for fields w/c^f ^ + Min(Af , A2) < he < 



fc^f 2 + Max(A^, A2). The spectrum for the low lying fermion, in contrast, shows all the 



desirable properties, vanishing at h^. = ^/A^A2 only for A; = 0; in addition, the effective 
gap, Ag ~ |/i — hc\, is as expected in the Majorana fermion representation of the critical 
Ising model and so is the relativistic dispersion for long wavelengths. Finally, when we 
integrate out the more massive fermionic modes we are left with a strictly non-interacting 
free Majorana fermion theory regardless of any zero-field interactions in ( |2.61 ); this is 
because all interactions will be polynomial in the one Majorana field left, and will vanish 
by fermi statistics. Thus we see that in contrast with the boson description, the free 
fermion theory is actually best near h = he- 

To summarize, on general grounds, one can expect qualitative agreement between 
both models up to magnetic fields close to he where the fermion model is expected to be 
a better description of the system. 

We now wish to look at some important matrix elements as phrased in the two models. 
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We start by defining, 



la,b{k,q) =< a,k\ l{0)\q,b > 



{21 



we use I = (j) X U for bosons to write 

wliere we define the cross product matrix with the Levi-Civita symbol by S* 



(2.89) 



— * — * — » — * 

Using, / = -^^(V'l ^ V'l + '^R ^ '4^r) fo^^ fermions, we write the analogous expression 



la,b{k,q) 



( 



(2.90) 



\ T^V^ SVR + R^aiV^ SVaiT T^V^ UV*T* + R^aiV^ i:V*aiR* 

X, ^, V,(7i,R and T were all defined in the sections on diagonalizing the models. As can 
be explicitly checked, the 0(3) free bosons are analogous to the NLcr model with the 
function, G{9), defined in Eqn. ( p.24|) , set to one. This is the general result for the 0{N) 
model for large A^, and makes sense, since the Landau- Ginsburg model is a large 
approximation to the NLo" model. In case of axial symmetry, one need only substitute 
the correct gaps into the energy factors: 



< k,a\r{0)\q,b>=ie" 



(2.91) 



with = ^A;2^2 + /^2_ 

The 0(3) fermion model exhibits a non-trivial G(^)-function. We can use the results 
from Eqns. ( glSBD and ( ^ ) in Eqn. {^M^ to calculate that 

1 



G{e) = -sech(^/2) 



{uk - k){ujq - g) + \/{uJk + k){uq + q) 



(2.92) 
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To obtain the U{1) results we, again, make the gap substitutions as done in Eqn. ( |2.91| ). 
This result is quite different than the boson prediction. It, in fact, vanishes with the gap 

— * 

for backscattering, k ~ —q. This is because the I does not couple left and right moving 
fermions while the opposite holds true with the bosons (and NLcr model). For small 
momentum exchange, all the models give universal predictions for matrix elements of 
S^{0). However, matrix elements of 5*^(0) at small momentum exchange are somewhat 
sensitive to the ratio of the gaps, in the boson model, while not at all so in the fermion 
model. In Chapter 6 we discuss experiments which might investigate this behaviour 
further. 

When the symmetry is orthorhombic there are few conservation laws to restrict the 
form of matrix elements of spin operators. Furthermore, when a magnetic field is added, 
Lorentz invariance is explicitly broken. We can, however, say that correlations among 
spin operators are still diagonal: This is true by virtue of the Z2 x Z2 symmetry. We also 
know that the new energy eigenstates, labeled by + and — , are mixtures of eigenstates 
of with eigenvalues = ±1. This guarantees that 

<k,+\S%->oc5i3 (2.93) 

It isn't terribly illuminating to write down the actual matrix elements. We can say, how- 
ever, that in the boson model, for h he, all matrix elements of form, < — , k\l^{0)\b, q >, 
which are not zero by arguments given above, diverge at /c = as fractional powers of 
{he — h). Everything is nice and finite with the fermions. This is another symptom of 
the sickness of the free boson model near criticality. Again we see that interactions are 
expected to play a crucial role in the boson description. 

We finish by describing some matrix elements near zero magnetic field. We expect 
that the intrabranch matrix elements, < ±, A;|S''^(0)|±, g >, vanish aX k = q with the 
field. For k = q = and h ^ 0, 
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<-,0|;»(0)|-,0>=-(^^-^j bosons 

2h 

— — fermions (2.94) 

where we've assumed, Ai > A2. The resuh for < +, 0|/^(0)|+, > is obtained by 
exchanging 2 and 1. Notice that this hmit does not commute with the U{1) hmit, 
A2 — > Ai. This is to be expected since these matrix elements are constant in the axially 
symmetric case. 



Chapter 3 



Model Predictions for Ti 



Let's recall the expression for 1/Ti, Eqn. ( |1.69| ), derived in Chapter 1 



Y = l^^ ^ -^f^i^ < {^i ^^'(0)} > (3.1) 

1 j,k,fi,u 

where h is taken to be in the z direction. As discussed in ( ^.2.21 ), only diagonal compo- 
nents of the spin correlation function will contribute. We also assume that the hyperfine 
tensor, Aij, is local 

-^ij ~ -A^ij (3.2) 

Thus we can write 

^ = Y1 l^^l' / dt e"'"^* < {S'{x = 0,t), S'iO)} > (3.3) 

where we've used translational invariance to evaluate the correlation function at the 
origin. We can now take a step back to Eqn. ( |1.68|) and write the above as 



1 ie-^-^'/T I p-^"/^! 
- = 27r 5: lA+'^n < ^'15^^(0)1^ > \H{En' - - ujnY- J i (3.4) 

^ n,n',u 

We will concern ourselves largely with the limit, un '^T <^ ^min (^ot^ that ljn ^ ImK 
for H ~ 1.5 T), so that the last factor in (|3^) can be set to 2e~-^"/'^. 

— * — * — * — » 

Consider now the operator in question, S{0) = 0(0) + 0(0) x iT(0). We wish to 
investigate whether dominant contributions to ( p.4|) come from the staggered field, </>(0), 
or the uniform part of the spin, 0(0) x i7(0). Let us first use the boson model to analyze 

53 
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Figure 3.1: First non- vanishing contribution to relaxation due to the staggered part of 
the spin 

the staggered contributions. This will be^ 

^ oc 5^ I < n'\<p\0)\n > \^6{E^: - E„ - a;;v)e-^"/^ (3.5) 

J 1 Stagg n,n' 

We assume, that we are in the regime T <^ ^min(^)' where Aj^jj^(/i) is the lowest 
(possibly field dependent) gap, or in other words, that the magnetic field is well below 
he, and we are therefore well justified in using the boson model (or NLa model) to describe 

— » 

the situation. Since 0(0) is a single magnon operator in the noninteracting theory, it 
only has matrix elements between states whose energies differ by a single magnon energy, 
Lj{k,h)] since is evaluated at the origin, k can be arbitrary. In particular, there are 
no matrix elements with energy difference, (which is essentially zero, compared to 
the other energy scales around). Including interactions, there will be contributions at 



finite T. The simplest process is shown in Fig. 3.1. It involves a -type interaction, as 



might occur in the Landau- Ginsburg or NLo" model. The vertical line represents the field 



"'^ There will not be contributions from cross terms between the staggered and uniform fields. These 
vanish because for < n|0|m >^ 0, one needs the number of magnons, n + m, to be odd, while this in 
turn implies < n\ l\m >= 0. Certain types of structural perturbations, such as discussed in Chapter 4, 
may change this analysis at high temperature and/or high fields. 
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(j). The incoming line from the right represents a thermally excited magnon of non-zero 
momentum, /c, and energy 2A (for simplicity, we use the isotropic model at zero field to 
make this argument; extending this to the anisotropic models is straight forward). The 
two outgoing lines to the left represent magnons at rest (recall that the wave vector, fc, 
is actually shifted by vr, and so the lowest energy antiferromagnetic spin excitations vary 
spatially as e^'^^). This diagram gives a non-zero matrix element proportional to A/A^ 
(where A parametrizes the 0^ interaction). Note, however, that since the initial and final 
state energies must be at least 2A, there will be a Boltzmann suppression factor of e~'^^l'^ 
to this contribution. Thus 

1 oc A^e-^/^ (3.6) 
^iStagg 

Including anisotropy and a finite field will give various contributions of this type. The 
greatest will be suppressed by exp(— 2Aj^jj^(/i)/T). It is also consistent to interpret this 
result as giving the single magnon a finite width at T 7^ 0. This, however, cannot change 
the conclusion that there is a double exponential suppression factor contrary to the model 



proposed by Fujiwara et. al. ||34] 



Let us now consider the contributions to 1/Ti from the uniform part of the spin: 

^ oc I < n'|r(0)|n > |25(E„, - E„ - c.^)e-^"/^ (3.7) 
J 1 Unif n,n' 

As discussed in Chapter 2, the 1-particle matrix elements selected above are non-zero 
in general, even in the non-interacting boson or fermion model. This is because /^(O) 
is a two magnon operator, able to create one magnon and annihilate another. In the 
presence of anisotropy and magnetic field, the three magnon branches are split, so we 



must distinguish between interbranch and intrabranch transitions (see Fig. |3.2|) . This 
is possible since contributions may come from all wave vectors. One set of important 
processes (ie. the ones corresponding to transitions between the lowest energy magnon 
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Magnon Dispersions For Lowest Two Branches 




Figure 3.2: Inter- vs Intrabranch transitions 



Chapter 3. Model Predictions for ^ 



57 



states) will come from single particle intrabranch transitions along the lowest mixed 
branch (intrabranch transitions are not allowed along the branch corresponding to = 
since h is parallel to z). This will be one of the leading effects with Boltzmann suppression 
of e"^"'-'*''/^, where A_(/i) is the lowest field dependent gap. It is important to realize 
that these will only be present if the hyperfine coupling, A"*"^, is non-zero. In fact the zero 
field gap structure and the choice of direction for placing the magnetic field will affect 
whether there are competing transitions. The next contribution, possibly as significant 
as the one just described, can come from intrabranch transitions along the second lowest 
branch and/or interbranch transitions between the lowest mixed branch and one of the 
other two branches. The Boltzmann suppression factor will, again, favour the lowest 
energy processes which occur at the gap to the highest branch involved in the transition. 
The important point is that the Boltzmann suppression factor for any of these processes 
is larger than that associated with contributions from the staggered field. To summarize, 
the dominant contribution to 1/Ti at T ^ ^min' ^^^^ come from Eqn. (|3.71 ). 

As we approach the critical field, the above analysis breaks down. As discussed 
previously, interactions are expected to become large in the boson model. Moreover, as 
the gap closes, the Boltzmann factors will fail to discriminate between the contributions 
of the uniform and staggered fields to 1/Ti. Arguments involving the fermion model 
are tricky because the staggered field has no simple representation in terms of fermionic 
operators. However, we explicitly show later in this chapter that the staggered component 
will dominate sufficiently close to hf.. 

Above the critical field, the analysis depends on the symmetry. For f/(l) or higher 
symmetry, the system remains critical and the staggered correlator remains dominant. 
For lower symmetry, the gap opens up once more, and sufficiently far above the critical 
field, we expect the uniform correlator to dominate once more. 
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3.1 Tf ' for h^hc 

In this section we concern ourselves with the regime discussed above, ojjq ^ T ^ 
^min(^)- calculate the relaxation rate in the isotropic, U{1) and Z2 scenarios. 
We begin by deriving a general result valid in this regime, and proceed to discuss its 
application in the different cases of symmetry. 

Consider a contribution to 1/Ti coming from transitions between branches r and s. 
Without loss of generality, we assume that r has a higher or equal gap to s (r, in fact, 
could be the same branch as s). We call the corresponding contribution to Eqn. ( |3.4| ), 
^ . This will be a sum over single particle states on r and s: 

1 =47r5:|A+f / ^5K(g)-.;,(A:)-^^)e~-('=)/^| < fc, r|r(0)|g, . > ^ (3.8) 
-'Irs i J (27r)^ 

Note that there will be a similar contribution with the labels s and r exchanged, if s 
and r are different branches. This corresponds to scattering an initial particle on the 
s branch through the hyperfine interaction to a final magnon on the r branch, or the 
reverse process. We take account of both of these possibilities later. Also, we are keeping 
lon finite to cut off infrared divergences which crop up in the intrabranch processes. We 
now do the integral over q to get 



rp / y I 



■X 



^'^^^KAk.Qim' + \uk,-Q{m') hrf (3.9) 

/ q=Q{k) 

where Q{k) is defined by uJs{Q{k)) = ujr{k) + un, and ^ is as defined in (|2.88|) . When 
Aj.(0) ^ T, the above integral will be strongly peaked at = 0. Moreover, we can 
neglect ujn- in ujr{k) +00^. The only factors in the integrand for which we should retain a 
^-dependence are the exponential and the possibly infrared divergent denominator, Q{k). 
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We therefore write 

-'-Irs i J q=Q(0) 



oo g-aJr(fc)/T 

It is not too difficult to show that to ffist order in small quantities, /c^ and a; at, one has 

^'i'^^-^'^o^^ml^J (3.n, 



We can also expand the exponent: 



ul,(k)IT : 

The relevant integral over k then becomes 



oo g 2tiJr(0)T 

dk , (3.13) 



fc2 



By changing variables, — ^^^^ > k'^, we can write ( p.l3| ) as 



dk , = (3.14) 



where 

«..(T,/.) = -^— ^^-^ (3.15) 

The /i-dependence of a^s will largely come from its dependence on C(Jr(0). The integral 
can be expressed in terms of special functions: 



e"-(^''^)/2iro(ar«(T,/i)/2) (3.16) 
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where Kq is the zero order modified Bessel Function. When the gap, uJr{0), is very large 
compared to the typical momentum, (5(0), exchanged in the transitions, (this is the case 
for intrabranch transitions), a^s — * 0. In this limit 

^aMT,h)/2j^^ («..(T, h)/2) ^ - log (a,,(T, h)/4) - 7 (3.17) 
where 7 = 0.577216... is Euler's constant. We can now summarize 

^ = E ^^^^^^^ (|/;,(o, Qmf + KM -Qmf) e--(°)/^ 

-'-Irs i ^ 



s 



e"-(^''^)/2i^o(«r.(T,/i)/2) (3.18) 



The full expression for the relaxation rate is 

^ = T.^ (3-19) 

1 r,s -'-Irs 

The effect of interchanging s and r in (p.8|) is therefore included in the above. 

An important thing to learn from the above calculation is that contributions from 
transitions between states involving small momentum exchange {Q —>■ 0) will dominate 
due to the logarithmic divergence in ( p. 181 ). This is particularly the case with intrabranch 



versus interbranch transitions. In intrabranch transitions one is allowed momentum ex- 



changes as small as Q ~ y2Ar{0)ujN/v. This will typically be much smaller than the 
smallest allowed interbranch momentum exchange, Q ~ (Ar(0) — As{0))/v. The conclu- 
sion is that, unless the branches in question are extremely close to each other, interbranch 
transitions will play a secondary role to intrabranch processes, even ignoring the more 
obvious suppression due to different Boltzmann factors. Of course, if the hyperfine inter- 
action has high symmetry, one will not see intrabranch transitions at all. This suggests 
that an NMR relaxation study could provide information as to the nature of the hyperfine 
tensor. 
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Before expounding on this result in the individual cases of different symmetry, we 
would like to mention the effects of higher temperature, or correspondingly, including 
the fc-dependence of the various terms approximated at = 0. In the more general 
Z2 X Z2 situation^, we are strictly justified in expanding the exponent in Eqn. ( p.l2| ) only 
for T <^ Ar(/i). For higher temperatures, one expects contributions from k > Ar{h)/v, 
where the expansion is not convergent; in this case one is better off numerically integrating 
( p.lOp (making the k = approximation for the other terms is still valid, as we will see). 
In either case, we can estimate the error in neglecting terms of order fc^". First, notice 



that all gaps are always greater than ^Jv'^k'^ + — h, where is the smallest zero 
field gap. We therefore write 

-ujr{h,k)/T 



■lk^ + Q^{k) J 



The last estimate is actually quite generous, especially for large n. In the worse case 
scenario of the Haldane phase, A/v ~ 1/4. This allows us to expect an error of at most 
10% in neglecting the fc-dependence of the terms in ( |3.10|) . 

We still have to estimate the error incurred in making the expansion in the exponential 
at T <^ Am. The next term in the expansion is — with = ^- This will 

give a contribution, 

J 8AIT ^fWTW 4A, ^ ' 

This is potentially more serious as T ^ A,, or A^ 0. To summarize, Eqn. ( p. 101 ) 
is generally a very good approximation; when T ^ A^, one can safely expand the 
exponential, while for T > A^, one is better off numerically integrating (|3.10|) . 



^One is more fortunate in the U{\) case; because of the siniphcity of the gaps, the expansion is good 
for k < A±, regardless of the value of h. 
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Figure 3.3: The gap structure for 0(3) symmetry. 
3.1.1 Isotropic Symmetry 

In the case of 0(3) symmetry, the field dependent gap structure is as in Fig. |373| . The 
lowest and highest branches correspond to magnons with = ±1 respectively. The 
middle branch corresponds to = 0. The interbranch gap is h. Using the result of 



the previous section and that of p.l.3| we can immediately write down the intrabranch 
contributions to 1/Ti: 

^ = lA-^r—.i^ogiAT/coj,) - 7](e-(^-^)/^ + e-(^+^)/^) (3.22) 
Ti Intra ttt;^ 

where in this case, Qq = 2ujnA/v'^. It is quite likely that higher dimensional effects may 
cut off this contribution at energy scales larger than uj^. For example, weak interchain 
couplings, Jj, would replace un in the above expression by a quantity of order J/. 
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The interbranch contributions between the lower two and higher two branches can be 
hkewise calculated to give 

1 = + |A+J^ni^e("+^)/2^/^o((/^ + -^)/2r)(e-^/^ + e-(^+'^)/^)(3.23) 

Ti Inter TTf ^ 2 A 



Following Fujiwara et. al. ||3J], we write 



1 = 1 +1 ^F(/.,T)e-(^-^)/^ (3.24) 
1 1 Intra ^ i Inter 

We see that the nature of F{h, T) depends largely on the form of the hyperfine coupling. A 

less general and somewhat more qualitative version of this formula was given by Jolicoeur 

and Golinelli [^], and by Troyer et. al. independently of our work. Jolicoeur and 



Golinelli discussed the isotropic NLcr model and derived only the leading exponential 
dependence on temperature; Troyer et. al. considered the Heisenberg ladder problem, 
which has a low energy one-magnon excitation spectrum identical to that in the isotropic 
NLcr model, and only included the leading interbranch transition in their expression. 

3.1.2 Axial Symmetry 

Here we are faced with two possible situations: the = branch can lie above or below 
the doublet. In the former case, the larger the interbranch gap between the doublet and 
the singlet branches, the more suppressed will be the interbranch contributions to l/Ti. 
On the other hand, in the latter scenario, inter- and intrabranch contributions will always 
be on the same footing (see Fig p.4|) . The expression for the intrabranch transitions will 
be essentially identical to the one in the 0(3) case: 

1 . .-.,24A 



■[log(4T/cu^) - 7](e-(^^-"^/^ + e-(^^+'^)/^) (3.25) 



Ti Intra ^rf ^ 

where A_l is the gap to the = ±1 branches and A3 is the gap to the = branch. 
The corresponding formula for the interbranch transitions is somewhat more subtle and 
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Figure 3.4: The gap structure for U{1) symmetry, (a) A3 > Aj_; (b) A3 < Aj_ 
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depends on the position of the singlet with respect to the doublet as well as on the size 
of the magnetic field. This is because ^ in (|3.18|) depends on which branch lies higher 
and the size of the gap between them. When the singlet sits higher, 

1 -(|A+^|2 + |A+^P ^ 



T^i Inter ' ' '" ' ' ttv"^ 

{e'^-^/^iro(«ir./2)e-^^/^A3(|/S_(0,Qi(0))p + |/-_(0,-gi(0))|^) 

+e"-»/2iro(«2.s/2)e-*^('^)/^M(/i) (|/S-(0, ^2(0))!^ + |/S-(0, -^2(0))^) } (3.26) 

where 

(A3 + - Aj Qijoy 

ttl rs = - 



2A3T 2A3T 

/i< A3 - A 



"2! 



(A3-fe)^-A^ ^ Ql{0)v^ 

2A3T 2A3r ^ ^3 



(Ax+fa)^-Ai _ Ql{0)v^ A ^ A - A 
2A I T ~ 2A I T ^ ^3 



A3 = M(h) /l< A3 - Ax 

M{h) = <( ' ^ ^ ' ^ (3.27) 

A± + h = M{h) + h h> A3- A± 

In the other case, when the singlet sits lower than the doublet, one has 

1 + ^ 



T^i Inter ' ' ' ' tcv'^ 

{e"-=/2iro(«i../2)e-(^^+'^)/^Ax(|/SL(0,gi(0))|2+|/o^L(0,-gi(0))|' 

+e'^-»/2iro(«2../2)e-*^W/^M(/i)(|/S_(0,g2(0))r + |/S_(0,-g2(0))r)} (3.2^ 
where this time, 

_ {A^ + hy-Ai _ Qijoy 

2A|T 2A|T 
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^2 



(Ax-fe)^-Ai _ QUO)-' < A , - A- 



(A3+fe)^-A^ ^ Q|(0)t;2 
2A3T SAsT 



h> A±- A. 



M{h) = { ^ ' (3.29) 

A3 = M{h) /i > Ax - A3 

A seeming catastrophe occurs when two of the branches cross at /i = IA3 — A_|_|. The 
interbranch contribution to 1/Ti diverges logarithmically. There are essentially two ef- 
fects that would cut off this divergence. Higher dimensional couplings can be counted 
on once more to replace Ql{0) as it approaches zero, with a quantity of order lOJj/J as 
derived in Eqn. ( [4.40|) in Chapter Four. Also, the divergence in the integrand leading 

to this problem is ~ —, i . This will be cured by a field with finite width. One 

still expects a peak in the relaxation rate, but this will be smoothed by the mentioned 
effects. 

3.1.3 Z2 X Z2 X Z2 Symmetry 

There isn't much more to say which would be illuminating in this case. We can, however, 

easily give the results for intrabranch contributions. These behave as the analogous 

expressions from the more symmetric situations. 

1 4|A+"i2 



Ti Intra 



-(log(4T/^^)-7)x 



|P_(0,0)p.;„(0)(^)^ ^e— + (--+) j (3.30) 
(0, 0) p depends on h as per Eqn. (|2.94|) . The formulae for interbranch transitions will. 



again, depend on the positions of the branches and the relative gaps between branches. 
Note that if there are intrabranch transitions allowed by the hyperfine coupling, then 
there will also be transitions between the + and — branches. Finally, from Eqn. (|2.94| ), 



we see that (13.30]) vanishes quadratically with the field. 
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3.2 Close to the Critical Field 

In this section we give qualitative results on the behaviour of the relaxation rate and 
in the process prove that the contribution from the staggered correlator becomes crucial 
as h he- We assume that we are now in the regime \hc — h\ T. In this limit, 
intrabranch processes along the lowest branch will dominate. Even if the hyperfine tensor 
possesses high symmetry (thereby ruling out intrabranch contributions from the uniform 
spin operator), we expect intrabranch contributions from the staggered part of the spin. 
Since the fermion model becomes exact in this limit, we will rely on its predictions. Long 
wavelength modes are now expected to play the most important role; we therefore write 
the dispersion relation of the lowest branch as 

ihc-h) + ^ 0(3) and Uil) cases 
Ivlk^ + A2 {Z2f case 



where vl = [aI+a% ' effective gap is Ag = 2{h — hc)hc/{Ai + A2). Now 

that the gap is actually smaller than the temperature, we must include multiparticle 
processes. This is simply done by replacing the Boltzmann weight by the appropriate 
occupation factors, ff{uj){l — ff{uj)). The derivation is straight forward and can be found 
in standard texts on many body physics (for example, see [§^). The uniform contribution 
to the relaxation rate is given by 

Due to the simple form of the density of states in the isotropic or axially symmetric case, 
one still obtains logarithmic behaviour for the above formula. In the anisotropic case. 



things are a bit different. We can combine the last expression with the results from p. 2. 2 
to get 

4|A+'=P A1A2 (u; + u;jv)sech^ 



(Tr%, ^ V. dk^Z^^I^^^iLMl (3.33) 
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At criticality we set Ag = 0. We may simply rescale the integration variable to obtain 

(Tr%^.fOcT (3.34) 

This is expected from the Ising model where the uniform part of the spin corresponds 
to the Ising energy density operator^, e, of scaling dimension 1. In terms of Majorana 
fermions this operator is ipL'4'R = ^- The correlator of the energy density operator on 
the infinite Euclidean plane is known from its scaling dimension and the restrictions of 



conformal field theory to be p4| 

< e(z)e(0) >= l/\z\^ (3.35) 

If periodic boundary conditions are placed in the time direction (corresponding to finite 
temperature), the correlator can be obtained by making a conformal transformation from 
the Euclidean plane into the cylinder (see -^p = ^^mzsT. 

< <M0) >= I ^^'^r'L (3.36) 

I Sm(i TTZ)\'^ 

Setting z = it + 6, we can get the contribution to Tf^ by integrating over / dte~^'^^^: 

fTfi),. .,oc [dte-'^^'--^j!pfi^^-— (3.37) 
V ^ /Unif J |sinh(T7r[t-i(5])|2 ^ ^ 

Changing variables, and assuming the integral is analytic as ujn (this can actually be 

proven by contour techniques), we see that by rescaling the time variable we reproduce 

Eqn. ( Pll ). 

We now turn our attention to the behaviour of the staggered correlator at criticality. 
We know the form of this function in both U{1) (and therefore 0(3)) and Ising cases 
from Ref. [Q. On the infinite Euclidean plane we have 

< 0t(;,)0(O) >^ |^|-| f/(l) 

_i (3.38) 

< a{z)a{0) >~ \z\ 4 Ising 



■^The rest of the analogy goes as follows: the magnetic field plays the role of temperature as is obvious 
from the form of the spectrum; the inverse temperature is analogous to the size of the system in the 
Euclidean time direction; finally, the staggered magnetization corresponds to the disorder field, a. 
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The field = 0^ + is the charged U{1) field of the boson model; there are no 
problems in using this (in the U{1) case) as long as we account for the interactions, a 
is the disorder operator of the Ising model. It is highly non-local in fermionic language, 
and aside from its dual, the order operator, and the energy density operator, is the only 
primary operator in the model. Once more, making a conformal transformation into the 
cylinder of circumference 1/T, 

< 0t(,)<^(O) >^ Uil) 

|sin(7r2T)|2 (3.39) 

< a(z)a(0) >~ Ising 

|sin(7rzr)|^ 

Once more, setting z — it + S and integrating over time, we get that in the experimentally 
important limit, T ^ a;jv, 

(^i"Ostag i^^TA/v')--^ + OM U{1) 
(^i"%tag ^ (^^/^«)"" + I^i^g 

(3.40) 

For both symmetries, this implies a significantly stronger contribution from the staggered 
part than the uniform part. In fact, as long as we are sufficiently close to the critical 
regime, perturbation theory tells us that the above result will only be suppressed by 
factors of order 0{\h — hc\/T). In order to observe this behaviour experimentally, one 
must have T sufficiently large (having a large anisotropy, Ai — A2, would also help), so 
that the decrease in relaxation with temperature is obvious. This would require that the 
experiment be done over a broad range of temperatures so that any constant contributions 
to 1/Ti could be subtracted. In any case, the above should at least serve to clarify that 
the staggered contribution becomes important in this regime. 

Farther still from criticality, the analysis breaks down. We do, however, expect the 
staggered spin contribution to infiuence Tf^ through to the region T ~ \h — hc\- 
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3.3 Above the Critical Field 

Far above the critical field, h — he ^ T, the situation becomes even simpler. In the 
0(3) and U{1) case the system remains critical. The relaxation rate will be dominated 
by the staggered part of the spin. The fact that develops a vacuum expectation value 
(or likewise, the non zero magnetization of the ground state) will have no effect on the 
relaxation rate since ujn isn't strictly zero: 

dt < 0\S\t,0)\0 >< 0|5^(0,0)|0 > e-^^'^* = 2n{M'Y6{uN)/L^ (3.41) 

where M is the magnetization. We can therefore say that for the 0(3) and U{1) models, 
the relaxation rate (assuming cuat ^ T) has the simple temperature dependence 

^ oc {2nA^T/v^r-' (3.42) 
-i 1 

where t] is the critical exponent of the staggered spin correlator. Haldane argued rj = 
I + 0(p), where p = \ M\/L [Q. Thus the field dependence of 1/Ti is only through rj. 

When axial symmetry is broken, the gap reappears for h > he- In this case, one can 
use the fermion model to calculate the relaxation rate. This is made much simpler since 
the gaps to the two upper branches are presumably much higher than the lower gap (by 
at least 2A_(/i)). Therefore only intrabranch processes along the lower branch need be 
considered. The result is 

- ' ' ^"^-/"^A- ?^ |e,_(0,0)P(log(4T/c.^)-7) (3.43) 



Ti TT V dk 



k=0 



At sufficiently large magnetic field, /i ^ Ai — A2, some of the expressions simplify: 

|/J^,_(0,0)|2 1 A_ ^ /i - (A2 + Ai)/2 
fdul\ 2v^ / A2 + A1 



The rate will drop exponentially with increasing magnetic field. 
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3.4 Summary 

We would like to summarize the main results of this section. At temperatures much lower 
than the lowest gap^, (this can be below the critical field or far above it in the case of 
Z2 symmetry), ~ e^^'^l'^ . This is due to the dominance of two magnon intrabranch 
relaxation processes. In the case of axial or higher symmetry, the only temperature and 
field dependence comes from 

^ (X [log(4T/c.^) - 7]e-(^-'^)/^ (3.45) 
^ 1 Intra 

Given such symmetry, this is a model independent result. Including anisotropics and 
interbranch processes, the relaxation rate is given by Eqns. ( |3.1(j| ) and ( ^.191 ). 



When the lowest gap is much smaller than the temperature, the dominant contribu- 
tions come from one magnon processes (due to the staggered part of the spin). When the 
field reaches its critical value, whereupon the gap vanishes, the relaxation rate is given 
by Eqn. (|3.40| ). Above the critical field, the system remains critical with axial symmetry. 



and the rate is then given by Eqn. ( p. 421) . With Ising symmetry, the gap reopens and 
eventually becomes large once more. 



^We refer to the lowest gap corresponding to a polarization direction perpendicular to the field. 



Chapter 4 

Material Properties and Possible Effects on Experiment 
4.1 Hyperfine Tensor 

Here we briefly discuss the effect of the nature of the Hyperfine tensor, Af^, on the NMR 
relaxation rate. As a reminder, /i and u are spin indices while i and j are spatial indices 
coupling the nuclear spin at site j to the atomic spin on site i. 

Assuming the magnetic field lies in the z direction, if A^j is isotropic in its spin 
indices, only |A^^p will contribute in Eqn. (|1.69|) . In particular, there will be no 
intrabranch contributions as these require a coupling to S^. Interbranch transitions 
will not be limited, but no contribution to them will come from the term proportional 
to A^^A . These statements still hold true for a hyperfine tensor diagonal in the 
Heisenberg spin basis. Note that this implies that for the 0(3) model, intrabranch 
transitions are prohibited so long as one assumes that the nuclear gyromagnetic tensor 
is simultaneously diagonalizable with the hyperfine tensor. 

In general, especially if the NMR nucleus does not coincide with the magnetic ion 
giving rise to the effective spin in the IDHAF, the anisotropics on the spin chain will 
not be simultaneously diagonalizable in the hyperfine tensor basis. Moreover, if there 
is more than one nuclear moment per spin contributing to the signal, it is unlikely that 
the effective Af^ will have the same symmetry as the nuclear Zeeman interaction. Thus 
conditions have to be quite convenient for intrabranch transitions to be missing from 
the rate. This can be important at very low temperatures where we can experimentally 
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distinguish between the processes. First of all, intrabranch transitions along the '— ' 
branch will increase exponentially with field (for example, in the 0(3) symmetric case, 
the behaviour is e~^^~'^^^'^). This can be a most obvious difference at low temperatures. 
However, as is clear from the discussion in the last chapter, the '— ' gap can lie above the 
= gap (for example, if one places the field along the a-chain direction in NENP); 
thus both inter- and intrabranch processes will feature the same exponential rise with 
field. Also, it is possible that the lower gap depends very weakly on the field for certain 
magnetic field directions. This is true for anisotropic materials where it is difficult to 
place the field along a direction of symmetry local to the chain. There are two ways 
to distinguish the transitions in these cases. The simplest solution is to repeat the 
experiment changing the magnetic field direction until one clearly sees the e~^~^'^^^'^ 
behaviour. Alternatively, one can try to extract information out of the low temperature 
behaviour, rather than field dependence. For T h, F(h,T) in Eqn. (|3.24|) will behave 
as log(T/co'7v) — 7 if intrabranch transitions are allowed. If they are prohibited, F{h,T) 



will more likely behave as ^ yT/{6 + h) where 5 is roughly the smallest interbranch gap 
at zero field. 

In principle, if one has enough information about the gap structure of the chain, it is 
possible to deduce the relative values of the hyperfine matrix elements from the relaxation 
rate. This may be done by comparing the ratios of the rates measured with magnetic 
field along each of the effective gap directions (or in case of high symmetry, any three 
perpendicular directions); assuming one has extracted the intrabranch contributions from 
the measurements of 1/Ti, one can then work backwards using 

^ oc (4.1) 
2i Intra 

(3 corresponds to the field direction) to arrive at ratios of the hyperfine matrix elements. 
This will be, presumably, model dependent even in the low field limit. In Chapter 6, we 
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suggest experiments which would distinguish between the models. 

Finally, we discuss nearest neighbour effects of Af^. When these exist, contributions 
to 1/Ti will come from the correlation function, 

< SmS^+^iO) >= ^E^-'™/^ < Sii{t)SUO) > (4.2) 

n 

where Sn is the nth fourier mode. Taking oo, we write 

< snt)suo) >=< smsno) > +/ ^ i^-'" - 1) < sms^.M > (4.3) 

Expanding in k, the odd contributions in k from the second term will vanish when 
integrated over all states. From the work done in the last chapter, we know that powers of 
/c^ will give small perturbations of order 2A^/f^. We see that including nearest neighbour 
contributions, the relaxation will be essentially given by making the substitution, 

in all previous expressions where, A^^, is the nearest neighbour hyperfine coupling. 
4.2 Impurities 

So far, we have dealt with a single, infinitely long spin chain. In real experiments, 
however, chains are always finite, and they come in three dimensional crystals, and so 
there are many chains of varying lengths in each sample. In this section we deal with the 
fact that these chains often end or have defects. This is what we mean by impurities. 
Of course one can introduce doping (for example, replace some Ni sites in NENP with 
Cu) to explore the issue further; we will restrict ourselves to 'pure' samples, although 
our treatment can be extended to doped samples. 

We start by describing a 'finite' chain. Exact work on a related S = 1 Hamiltonian 
- the 'valence bond solid' [^], which also features a unique ground state (in the ther- 



modynamic limit) and a gap - has indicated that at the ends of a finite chain there are 
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free spin-| degrees of freedom. It is expected that the same holds true for finite chains 
of the Heisenberg antiferromagnet. Indeed, there has been convincing numerical and 



experimental evidence for this conjecture The same evidence also supports the 



notion that these nearly free end spins can interact with each other, essentially exchang- 
ing virtual NLo" bosons. The interaction is exponentially decaying with the size of the 
chain: 

Hj oc (-l)^e-^/« So ■ Sl (4.5) 

Here, the correlation length, ^, is roughly six lattice spaces. For experimentally realizable 
lengths, this interaction is negligible and the finite chain has a fourfold degenerate ground 
state. 

In a 'pure' sample, these finite chains will lie end to end, or be separated by some 
non- magnetic defect. It is therefore reasonable to presume that two adjacent end spins 
will interact with an effective exchange coupling, J', that will vary in strength from zero 
to something of order J. Recent work has explored this situation in depth. For weak 
coupling between adjacent end spins the effective Hamiltonian can be written. 

He = J' Si - §2 = a^J' s[ ■ S^ (4.6) 

-►/ , 

where S is a spin-^ operator and a is the projection of a spin-1 end spin into a spin- 
^ subspace. From numerical work |41, 16], we know that a ~ 1. This immediately 



gives a low lying triplet above a singlet with a gap AE = a'^J' (we assume that J' is 
antiferromagnetic. The ferromagnetic case is expected to give similar results, reversing 
the order of the singlet and triplet, but numerical work has not yet been done to support 
the analysis in this limit). This triplet will sit inside the Haldane gap. The triplet 
corresponds to bound states at the chain ends; this has been seen numerically in 



which also demonstrated that the above first order perturbation theory result for AE 
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_ tl i— |ti>-|lt> 

Figure 4.1: Impurity level diagram when D' = 0. 

is accurate up to Je ~ -37. The embedded states become delocalized and join the 
continuum at about J' ~ .7 J. This picture is unchanged up to about J' ~ 2J, after 
which the triplet returns as a bound state in the Haldane gap to merge with the singlet 
state as J' —>■ oo. In the type of chain we are considering, J' is unlikely to become much 
greater than J; it is much more likely that defects in a pure sample will serve to reduce 
the effective coupling between sites rather than enhance it. 

To understand the effect on NMR relaxation we explore the environment of NMR 
nucleii near the end spins . Take, for example, the nuclear spin coupled to S^. It sees the 
Zeeman split level diagram shown in Fig. [4.1| . If we assume that the 'free end spins' are 
not completely free, but are weakly coupled to the magnons on the chain, then relaxation 
can occur in two ways: when two levels are Un apart the chain end spin could make a 
transition - this will happen for a fields /i ~ J'. Alternatively, a thermal magnon coupled 
to the chain end could decay into another magnon with or without the accompaniment 
of an end spin transition - again, the energy difference between initial and final states 

^We would like to thank D. MacLaughlin for privately communicating his suggestions on the effects 
of end spin excitations in NMR. 
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must be t^AT. Marked on the diagram are the transitions that can be induced by S : sohd 
arrows represent possible transitions that potentially do not require coupling to the rest 
of the chain; dashed arrows represent transitions that could occur only if the impurity 
is coupled to the rest of the chain; solid circles represent spectator transitions which, 
again, require magnon assistance. It is easy to see that in all the above scenarios, the 
transition will be broadened by thermal magnons. This implies that the characteristic 
width will have a temperature dependence exponential in the Haldane gap. There is an 
additional mechanism, which we now discuss, that can, in general, affect this picture. 
If the material is anisotropic, with for example, a D-tjpe anisotropy, we must add the 
following term to the chain-end effective Hamiltonian, Eqn. 



He^He + a^D ((^f )2 + {S^')^) + D'SfS'^ (4.7) 

The first term is the familiar on-site anisotropy; it will only contribute a c-number to the 
effective Hamiltonian. The second term is allowed by symmetry, and we presume that it 
is a consequence of the defect (which arguably, would manifest itself in accordance with 
the available symmetry). We assume D' > — J' , so that the exchange interaction is still 
antiferromagnetic in the z-direction. As a result of this, the two = levels will shift 
by — while the = ±1 levels will be shifted by The transitions induced by the 
hyperfine interaction are shown in the new level diagrams in Fig. ^.2| . It is this more 
general case for possible transitions which we now carefully analyze (further anisotropic 
purturbations will not qualitatively change the picture). 

We begin by characterizing the interaction between the bulk magnons and the end 
spin. Since we have no information as to the nature of this coupling we will parametrize 



it in the spirit of Mitra et. al. using free bosons which after scattering with the 
impurity obtain a phase shift. We make two assumptions: first, that leakage across the 
impurity site is negligible, and second, that the impurity spin coupling to the boson does 
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ltt> 



\u>- ln> 



Figure 4.2: Impurity level diagrams for D' ^ 0. 

not allow for exchange of spin (this is consistent so long as J' is sufficiently weak) ; bosons 
on each side of the impurity will have wave-functions of form: 



(4.8) 



where i refers to the boson branch, k is assumed positive, |C±(A;)p = 1/2L and ± refers 
to the sector of the impurity spin with s^' — ±|. The boundary condition, 0^j(iv) = 0, 
gives 



= {nn - 5l{k))/L 



(4.9) 



We assume that the phase shifts are small and grow with increasing energy. This is 
tantamount to assuming a large step potential barrier of infinite extent (thus allowing 
no leakage for states below the barrier). This is certainly true in the limit J' — > 0. A 
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heuristic ansatz which has this behaviour is 

6i{k) ~ (4.10) 

Notice that the 0"*"'s need not be orthogonal to the ^'"'s, since these states are in separate 
Hilbert Spaces. The energy of one of these bosons at low temperatures is given by the 
free form 

Ei:^^A, + ^^ (4.11) 

Before going on we note that a;jv is typically much smaller than the energy level spacings 
due to finite size, for typical chain lengths. For example, in NENP, SE ~ .04 Kelvin for 
chains L ~ 1000a. One may therefore question the validity of boson assisted transitions 
when the bosons lack the ability to 'fine tune' a transition so that the difference between 
initial and final states is a; at. What saves the day, in this case, are higher dimensional 
effects. For sufficiently long chains or sufficiently strong interchain couplings, these will 
densely fill the spacings in energy levels along the chain direction. Assuming this is the 
case, (we show the conditions for this explicitly in the next section) we will not worry 
about this point further. 

Consider the coupling of a nuclear spin to one of the impurity spins, say S[ (from 
here on we will implicitly write S[ — S', for ease of notation). The familiar formula for 
the transition rate can be cast as 

1 ^ T J ~oo Till 

^ ^ M n,l,n',l' 

+ < ni;l\e-^^''''+''^^S'^'\n'ir,l' >< n[r,l'\e''^''''+''^^ S^'e-''^''''+''^^ni;l >} (4.12) 

where /j, = —jl = ±, 0. Hf, is the free boson Hamiltonian; |/; mi > denotes a state of the 
impurity spin (ie. an eigenstate of He), which for brevity, we denote as |Z >; n; denotes 
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the boson content of the state; in general, |n;; / > will contain two different multiparticle 
free boson states that have projections onto the Sf — ±| subspaces of \l;mi > (with 
appropriate phase shifts). To elaborate and make things a bit clearer, take the = 
state from the triplet. 

|ni,o; 1, (l T, i> ®|n+ > +1 i, T> ®|n- >) (4.13) 

where the states \n^ > correspond to the phase shifted bosons with Sf = ±|, respec- 
tively. Making the approximation, o; = 1, we return to the relaxation rate: 

1 2 

— = - / (it e-''"~*^|A-'^|2 

Til -2 ^-oo ^ 

E e-^(^"'+^')cos(tK, - E^,^+ei - e,])\ < l\S'''\l' > H < ni\n[, > \^ (4.14) 

n,l,n',l' 

Let's pick a particular transition and work it through. Consider /x = +, |/ >= |1,1 > 
and \l' >= |1, >. This transition could be of the type we've been discussing where for 
h ~ D'/2, we expect strong resonance if D' > 0. The expression for the rate becomes 

1 (z,r,/i) = 4e-''(-^'^^'/2-/.) r dt e-'^'-'iA-^i' 

J-IJ 2j J — OO 

Re |e^*(^'/2-?.) ^ g-/3£;„+g^t(£;„+-s„_,)| ^ ^+|^-, ^ J ^^^^^^^ 

I, n.n' ) 

Since the boson multiparticle states, < n|, are direct products of symmetrized free N- 
particle states, and since the energy of such a state is the sum of single particle energies, 
we can write the last equation in terms of single particle states: 



-OO 



Re|e^W2-ME^p/'_2^+ ^ e-^K+^'K+-K-'\<n+-i\n-'-j>A\ (4.16) 
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where we've used the fact that the grand partition function for noninteracting bosons 
is the exponentiated partition function for a single boson, n and n' now index single 
boson states, and i and j denote boson branches. The overlap of the boson states can be 
calculated from the form given by Eqn. (|4.8|) , 

(4.17) 



+ ., / . ,2 1 e /sin(5i(A;) + 5i(A;')) siMik) - 5Uk')y 



where we have parametrized the momenta of n and n' with k and fc', respectively. Finally, 
we write the one particle partition function (we can ignore the phase shifts for this 
purpose) as 



n,i 



-I3E'_ 



< n^; i\n :i > P 



Combining all of this allows us to write the exponential in Eqn. ( 4.16|) as 



-Ai{h)/T I dlidk^^ 



k'k 



i \2 



A;2 - k'^ 



(4.18) 



xe 



-1) 



(4.19) 



E- 



-^>W/^^(AV-AL) 



oo 







J{x-x')Tt 



= 0(t) (4.20) 



{x — x'Y 

The effect of the phase shifts is contained in the factor, (A^ — A*!)^, as seen from Eqn. 
( [4.10| ). Corrections to this due to Oijc") contributions to 5\.{k) will be suppressed by 
factors of 2AjT/f^. The imaginary part of 0(t) will shift the resonance from h = D'/2. 
This shift, at low temperatures, will be negligible. We are interested in the long time 
behaviour of 0(t). In this limit, the real part of 0(t) becomes: 



Ree(t) = - ^ e-^'('^)/^^^(AV - Xlf dx dx'e 



.si\i^{{x~ x')Tt/2) 



[X — X' 



i\2 
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The expression for the relaxation rate becomes 



X / rfte-*"^*cos(t(D72-/i))e-l*l^(^) (4.22) 



(4.23) 



I ^-/3{j'+D'/2-h) r(r) 



r2(r) + (/i-z}72)2 

This is the most important equation of this section. The other transitions can be treated 
the same way to arrive at analogous results. The key issue to note is that the impurity 
relaxation rate is an extremely sensitive function of the temperature and field. At tem- 
peratures well below the gap it is essentially a delta-function of h. As the temperature 
increases and becomes comparable to the gap, the rate broadens rapidly. 

Before summing up, we discuss the other possible transitions. First, notice that 
changing /j, has the same effect as reversing the sign of h and exchanging / and /': 

II fi < — > h ^ —h < — > I ^ I' (4.24) 

Note that in this simple model of boson-impurity coupling there are no transitions via 
S^', and therefore no transitions between the singlet and the = state of the triplet. 
In other words, the solid circles in Figs. [4.1| and ^4.2| are ignorable as are the dashed lines 



from |1, > to |0 >. This is expected in all but the most extreme of anisotropic exchange 
impurity models. Furthermore, the effect of reversing the spin states on the triplet is the 
same as reversing the sign of magnetic field: 

mi, mil — > —mi,— mi' < — > h —h (4.25) 
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Finally, the result of exchanging the |0, 1 > state with the singlet amounts to adding J' 
to the associated energy factor in the Lorentzian. A final expression for the impurity 
relaxation rate involving all eight possible impurity level transitions is 

^11 i=l f=l i U J + ^fi 

where Ei denotes one of the four possible initial impurity states, and Eji is the difference 
in energy between the initial state and one of the two possible consequent final states. 
The factor of two represents the contribution of both end spins on each chain (we neglect 
surface effects) . Now it can be seen more clearly that all but two of the elements in the 
sum above will contribute little due to the narrow gaussian form. The important terms 
are those where the energy in the gaussian is small; this can happen for certain magnetic 
fields: h ~ |-D'/2| and h J' + D'/2. In Heisenberg chains we might expect J' ^ \D'\ 
for most defects. Furthermore, the impurity contribution should be most evident at 
lower fields where the gap still lies high. Consequently, in experiment, one expects the 
h ~ |-D'/2| transition to dominate the picture of impurity contributions to 1/Ti. 

In a real sample, the NMR signal from the impurity will be proportional to the density 
of the impurities. Moreover, since defects will vary from chain to chain, one would be wise 
to average over a random distribution of couplings, J' and D'. In practice, experimental 
data could be analyzed for the 'peak' values of J' and D'. One could then model the 
distribution of couplings with the appropriate peak values. This could, in principle be 
checked against low temperature ESR measurements which ought to concur with the 
impurity model. 

A final expression for the relaxation rate due to impurities is 

(^) ^nj dJ' dD' p{J' - J')p{D' - D')^ {f, D') (4.27) 
V 1 / Imp •' J- 1 1 

where n is the density of impurities (or inverse length of the average chain); p is some 
distribution function. 
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4.3 Interchain Couplings 

In previous sections we mentioned tlie effects of intercliain couplings on various aspects 
of the physics. We now examine these in more detail. 

Nearest neighbour interchain couplings will enter the Hamiltonian as 

H ^ H + Jj J2 (4-28) 

<i,j> 

where < i,j > index nearest neighbour spins not on the same chain. We can return to the 
derivation of the NLa model to see the effect of this additional term. Taylor expanding 

— * 

the continuum representation for Sj and assuming reflection symmetry about a site, 
Eqn. ( |1.38D will change to 

where we chose the z-direction to be along the chain, and the vector, <^j, is the dis- 
placement vector to the ith nearest neighbour of a spin not on the same chain (again, 
we assume that | CI is smaller than the correlation length. Note that the correction to 
dynamical part of the Lagrangian will correspond to J — > J + J/. Presumably, Jj <^ J, 
meaning that we were justified in ignoring this term. Setting the lattice spaces. Ax, Ay, 
to 1, we can now write an effective Landau- Ginsburg Hamiltonian to describe the physics; 
Eqn. ( p.. 50 ) will read 



The leading relevant interaction terms will always be local. Ignoring these, the resulting 
equations of motion are 

= [vdl - — + 2J:s Q ■ V)^) $ (4.31) 
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The dispersion relation becomes, 

uj^ = v^k' + A^ + vl C^ ■ kf (4.32) 

i 

where v±_ oc ^fJjJ . It is from this last formula which we now extract qualitative infor- 
mation about interchain coupling effects. 

First, recall that we claimed that for finite chains of certain lengths in real experi- 
mental situations, we no longer need to concern ourselves with 1-D finite size effects. In 
other words, we said that energy levels arising from interchain couplings will densely fill 

2 2 

the small gaps between magnon energy levels, AE" ~ -^^12 ■ Let's calculate this length in 
terms of J and Jj. We start by assuming a simple form for the interchain contribution 
to the dispersion: 

^iE(C- kf^vlal{kl + kl) (4.33) 

i 

where is some typical interchain distance, and expected to be 0(1). The size of the 

2 2 2 

interchain band will be ^-^^^ . Setting this equal to the gap in the magnon levels we get 



J 

,2 



~ ^ = 5- (4.34) 



In NENP, for example, this corresponds to lengths of approximately 100 lattice units. 

There is also the issue of cutting off divergent integrals which we discussed in Chapter 
3. In calculating transition rates, one often encounters integrals such as 

j dkdq5{uji-uji-E)f{k) (4.35) 

When E is close to the gap between the branches, uj]. and o;^, this integral can diverge 
logarithmically in the infrared. If one introduces interchain couplings, the integral over 
the delta function becomes 



(27r) 



Kzjdq d^kx. d^q±S{ui -uj{-E) 
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^ (47r)2ai ^ 2Ai 2Ai 2Aj 2Aj ' ^ ' 

where we have assumed a simple form for the interchain dispersion. For ease of calcula- 
tion, we now assume that Aj = Aj = A. The integral becomes 

dk_^ dq^ /)/(2i 2(2 2 2\ 



v'j_kj_ — v'j_q'j_ 

A 



J dkl^v^k^ + vlkl (4.37) 



At low momentum, k, where we need a cutoff, this integral is approximately 

Att 



6v±va± 



If we write the integral in Eqn. ([4.35| ) as 



(4.38) 



'4fdk^^ (4.39) 

then the cutoff, C, is seen to be 

c . lilHM ^ ,,J, „ 

We recall that the for intrabranch transitions, Qo ~ 2AuJi^/v^. Comparing this to C, we 
find that the cutoff becomes important for ujn < 70 J/a^. For example, in NENP, where 
Ji ~ 25 mK, we expect the cutoff to significantly dominate over the Larmour frequency. 

4.4 Crystal Structure 



When analyzing experimental data in terms of the idealized Heisenberg model with on- 
site anisotropies, one must keep in mind that the symmetry of the proposed spin-chain 
Hamiltonian may be constrained by the symmetry of the crystal and the local symmetry 
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about the magnetic ion. Additional terms may be added or subtracted to accommodate 
the structure of the substance, and these can have a great effect on the interpretation 
of data. Some important questions which must be asked before deciding on a model 
Hamiltonian for the material are: is the local crystal field symmetry about the magnetic 
ion commensurate with the symmetry of the unit cell? Is there more than one chain 
per unit cell? If so, are all chains identical? Is there more than one magnetic ion of a 
single chain per unit cell? If so, is there true translational symmetry from one spin site 
to another? 

In the next chapter, we analyze experiments performed on NENP. In so doing, we 
will address such considerations. 



Chapter 5 



NENP: Direct Comparison with Experiment 



5.1 The Structure of NENP and Experimental Ramifications 



A schematic diagram of Ni{C2HfiN2)2N02{ClOi) (NENP) is shown in Fig lOj. Each 
chain is comprised of Ethylenediamine-Nickel chelates separated by nitrite groups. The 
magnetic ion is Ni^^\ experiments indicate that these ions interact antiferromagnetically 
along the chain with coupling J ~ 55K. There is a large single ion anisotropy, D ~ 12K, 
as well as a small axial symmetry breaking anisotropy E ~ 2K. Interchain couplings are 
estimated at J// J ^ 10^'' 

It is important to realize that two neighbouring Ni"^^ along the 6-direction are not 
equivalent; rather, one is related to the other by a tt rotation about the h axis. Also, the 
angle along the — Ni — O bond is not exactly vr, meaning that the Ni site is not truly 
centrosymmetric. Most importantly, the local symmetry axes of each Ni ion are rotated 
with respect to the ahc (crystallographic) axes. To demonstrate this we now note the 
coordinates of the Nitrogen atoms in the Ethylenediamine chelate surrounding the Nickel 
(placing the Nickel at the origin): |^ 



Atom 


a (A) 


6(1) 


c(A) 


N{1) 


2.053 (3) 


.162 (3) 


.338 (3) 


N{2) 


.619 (3) 


-.184 (3) 


-1.971 (3) 
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Figure 5.1: NENP 
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Figure 5.2: Local and crystallographic axes projected onto the ac-plane in NENP 



The other Nitrogen atoms in the chelate can be obtained by reflection through the 
Nickel. One easily sees that projecting this structure onto the b plane yields symmetry 
axes (in the b plane) rotated ~ 60° from the ac-axes. This is shown in Fig. |5.2| . The 
inclination of the local Nickel axes from the abc system can be obtained by taking the 
cross product of the two Nitrogen vectors (ie. the normal to the plane described by the 
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four Nitrogen atoms in the chelate: 

ri = (-.06 (1), .98 (1), -.11 (1)) (5.1) 

The local Ni 6'-axis makes a ~ 10° angle with the 6-axis, while the azimuthal angle 
in the ac plane is ~ —28° from c. The 10° tilt is roughly about the a' direction of the 
local symmetry axes. 

One may worry that the NO2 group may distort the local symmetry axes, but re- 
markably enough, when projected onto the ac plane, the three atoms in the molecular 
ion sit on the c' axis. This reinforces our suspicion that the local symmetry axes are 
indeed the above. 

Next, we consider the whole space group of NENP. The most recent attempt to solve 
for the crystal symmetries has concluded that the true space group of the material is 



Pn2ia [45|; this is a non-centrosymmetric space group with a screw 2i symmetry about 
the b axis, diagonal glide plane reflection symmetry along the a axis, and an axial glide 
plane reflection symmetry along c. Experimentally, attempts to solve the structure in 
Pn2ia have not been successful; rather, it seems that Pnma gives a better fit. The 
main difference between the two is the presence in Pnma of a mirror plane parallel to 
b at centers of symmetry at various locations in the unit cell, and two-fold screw 
axes separating these centers of symmetry. The reason for the experimental discrepancy 
is attributed to disorder in the orientation of the nitrite group, the perchlorate anions, 
and the existence of a local or pseudo center of symmetry lying very close to the Ni 
(thousandths of an Angstrom) A crucial point is that both space groups share the 
axial glide planes along a, the diagonal glide planes along c and the 2i screw symmetry 
about b. These generate a total of 4 Ni sites per primitive cell and two chains through 
each cell. The two chains are such that the Ni chelates on one are the mirror image of the 



other. Figure shows a projection of this picture onto the ac plane. The presence of the 
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Figure 5.3: A projection of the NENP unit cell onto the ac-plane showing two chains per 
unit cell 
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2i screw symmetry about each chain axis introduces staggered contributions to the local 
anisotropy and gyromagnetic tensors. This is because, as motivated above, these are not 
diagonal in the crystallographic coordinate system. The resulting spin Hamiltonian is 

H = JJ2[Sr + ■ T>S, - fiBS. ■ GH + (-1)^(5, ■ dS, - /i^S, ■ gH)] (5.2) 

i 

We make the assumption that the symmetry of the anisotropy and (yf-tensors is the 
same (ie. that at each site they can be simultaneously diagonalized) . This is rigorously 
true when the crystal field symmetry about the magnetic ion is no lower than orthorhom- 
bic (a sketch of a proof is found on p. 750 of |Q ). We can get the required parametrization 
for the (yf-tensors from high temperature uniform susceptibility data . This is based on 
the idea that at high temperatures the Ni atoms will behave as an ensemble of uncoupled 
spins (s = 1) with the same gyromagnetic tensor as in the antiferromagnetic case. With 
this in mind we get 



G 



^ Gc'Cos\e) + Gb>sm\e) 





V 





Ga' 

Gb' cos\e) + Gc' sm\e) ] 

\ 



(5.3) 



sm{e)cos{e){Gb' -Ge, 

(5.4) 

\sm{9)cos(9){Gb' -Gc') j 

Here 6 ~ 10°, and Ga' = 2.24, Gb' = 2.15, G^' = 2.20 are the values for the local 

G-tensor that give the observed high temperature gf-tensor when averaged over the unit 

cell. Correspondingly, the anisotropy tensors must have the following form: 

/ n n n \ 



D 



D^, 

Da' 

Db' / 



(5.5) 
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tan(29) , 



(5.6) 



V J 

The parameters Da>,Dc',Db> are to be fitted by experiment to the model used to 
describe the system. The boson Hamiltonian can now be written 

H = J dx[^ n + + ^ H -Gi^x n)+ (5.7) 

i ■ d( X n)-fiB$-gH + X{ ^V] 

V 

The term containing d breaks the Z2 symmetry along the a' (lowest mass) direction. 
It will also renormalize the masses. The second effect can be ignored in the approximation 
that the (p^ term is ignored if we assume the masses are physical. Symmetry breaking, 
however, leads to the presence of a static staggered field even below a critical magnetic 
field. A gap will always persist. The staggered field term will break the Z2 symmetry 
along the c' or b axis, depending on whether the field is applied in the b or c' direction, 
respectively. A static staggered moment will likewise appear due to this term. The 
effect on the relaxation rate will be small, although there may be consequences in other 
experiments [^, ^ . 

We would now like to discuss the effect of having two inequivalent chains per unit cell, 
with local axes different from the crystallographic axes. We label the two chains found 
in a unit cell of NENP 'chain 1' and 'chain 2' corresponding to the chains in the upper 
left and lower right corners of Figure |5.3| respectively. The dispersion branches of chain 1 
are given by Eqn. (20) of (the expressions are roots of a complicated cubic equation 
and we feel that citing them will not prove illuminating) only the field is a — 30° from 
the c' axis where a is the angle of the field from the crystallographic c-axis. Similarly, 
the dispersion branches of chain 2 are calculated with the field a — 150° from the c' axis. 
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Experiments which average over signals, hke susceptibihty or NMR Ti measure- 
ments, must consider their resuhs an average of two different measurements (correspond- 
ing to the two different chains with their relatively different applied field configurations). 
On the other hand, experiments such as ESR, should show a separate signal for each 
chain. The NMR relaxation calculations performed in Chapter 3 assume the field is 
placed along one of the crystal axes. In this special case, the dispersions for the two 
different chains are identical. Although the dispersions will be more complicated as will 
be the matrix elements, l^hiO, 0), we do not expect great qualitative differences between a 
calculation as done in Chapter 3 and one which accounts for the actual symmetry when 

— » — * 

the field is placed along a crystal axis. There will also be contributions due to the (p- I 
correlator; these are also expected to be small. There are, however, important manifes- 
tations of having two inequivalent chains. These will be discussed in the next chapter 
when we suggest further experiments. 

In conclusion, consideration of the crystal structure introduces both symmetry break- 
ing terms and two inequivalent chains per unit cell. The symmetry breaking terms will 
give small corrections to the relaxation rate. 



5.2 Analysis of the Data 

By far, the most studied Haldane gap S = 1 material is NENP. The most recent measure- 
ments of the relaxation rate, 1/Ti, on this substance have been made by Fujiwara et. al. 



3J] . Before directly comparing our results to the data we discuss the expected results on 
a pure (infinite) system. The three gaps are given by neutron scattering: = 1.17meV, 
Af, = 2.52meV and = 1.34meV. We use v = 10.9meV, and the generic value of 2.2 for 
the electronic gf-factor. Since we do not have an accurate description for the hyperfine 
coupling of the Ni ion to the protons in its surrounding chelate, we assume a uniform 



Chapter 5. NENP: Direct Comparison with Experiment 



96 



value for all the contributing hyperfine matrix elements in a given direction of the applied 
field. Writing 

Tf^ = F(/i,T)e-^-('')/^ (5.8) 
we use the results from Chapter 3 to plot F(h,T) for bosons and fermions and for fields 



along the chain a, b and c directions. The results are shown in Figs. |5.4| - pTi . 

We included multiparticle transitions by simply replacing the Boltzmann factor by 
appropriate occupation factors in Eqn. ( |3.18| ): fb{^ + fb) = cosech^(^)/4 for bosons, and 
//(I — //) = sech^(^)/4 for fermions [Q. Within approximations used, multiparticle 
effects amount to multiplying the final expressions by (1 ± e'^'^=)~^. At higher tempera- 
tures it is also necessary to include the fc-dependence of the integrand past the peak at 
the origin. We expect that at temperatures T ^ ^ fields /i ~ ^ the numerically 
integrated results would differ by about 10 percent. 

F{h,T) is shown for fields up to 9 Tesla even though the {(3uj- ^ 1) approximation 
is no longer valid at such fields. This is done to contrast the predictions of the boson 
and fermion models. It's easy to see that the boson result for F{h,T) diverges at the 
critical field, while no such catastrophe is present in the fermion result. This divergence 
is logarithmic and infrared. It will persist even after account is made for the staggered 
part of the correlation function. Multiparticle scattering will in fact worsen the effect, 
since the bose distribution function diverges as l/cu with vanishing energy u. This again 
is evidence of the inadequacy of the free boson model close to criticality. 

In NENP, when the field is along the b direction, we expect relevant interbranch tran- 
sitions only for small field. In this regime, one must also be careful to include intrabranch 
transitions in the second lowest branch. All these processes are of the same order. Even 
though the intrabranch rates vanish at low fields, the interbranch contributions are sup- 
pressed by the absence of low momentum transitions (ie. Q for the interbranch transitions 
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is 0(yAf — A2) as opposed to 0{uj]y).) For this case, only need be calculated. 

When the field is along the c direction (corresponding to the middle gap), we restrict 
ourselves to calculating intrabranch transitions along the lower branch and interbranch 
ones between the lower and c branch. There are no intrabranch processes along the c 
axis. Calculating the interbranch transitions amounts to calculating l/^cP l^acP- 
When the field is along the a direction, the calculation proceeds as above. The crossing 
of the branches provides for the interesting effect mentioned earlier. The peak in 1/Ti 
can be used to locate the true ac-axes for the chain 

Notice that F{h, T) for the field parallel to the b axis is nearly field independent over 
a large range of the magnetic field. This behaviour is quite easy to understand from 
the universal results valid in the axially symmetric case, discussed in Chapter 3. When 
the field is along b, the system is only slightly anisotropic, and so the axially symmetric 
results roughly apply. Fb is roughly independent of field with axial symmetry since 
is nearly h independent (in fact, Ff, exhibits a logarithmic divergence as /i — 0). On the 
other hand, Fc vanishes quadratically as /i — > 0. Including the small breaking of the axial 
symmetry corresponding to Ac — = 2°K, Ft is essentially constant down to low fields 
of order A^ — A^ ^ IT, before rapidly decreasing as seen in in the figures. 

We now proceed to directly compare our results with those of Fujiwara et. al. Since 
the hyperfine coupling is not known, we find a best fit to it using the experimental data. 
This is best done for mid-sized fields: in the low field regime impurities may dominate, 
and in the high field regime the staggered part of the spin is expected to contribute. Figs. 
|5.8| and ^.9| are such fits to the boson and fermion models. 

In producing these fits we get different values for Ay^, the hyperfine coupling for a 
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Figure 5.8: Theoretical (lines) vs. experimental data (circles and squares) 
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Figure 5.9: Theoretical (lines) vs. experimental data (circles and squares) 
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field placed along the i(7-direction: 



8.5 MHz fermions 
Ab ~ { (5.9) 

7.0 MB.Z bosons 



17.7 MHz fermions 
A, ~ { (5.10) 

12.0 MH^ bosons 

These values are reasonable for dipolar hyperfine couplings between a nuclear spin (s = 
1/2) and the Ni spin at a distance of about 2 A: 

yl^/^^2MH^ (5.11) 

Also, we can get a similar feeling for the size of the hyperfine couplings from Knight shift 
| 51| and magnetic susceptibility data for a field placed along the 6-axis. 



IXnSH ~ AfXbXsH (5.12) 

At about 4K, the susceptibility is roughly a fourteenth of its maximum value. Given that 
Xyiax ~ that the Knight shift at large fields is about 10"'^, we get A ~ 8 MH^;. 

It should be noted, however, that these are order of magnitude estimates; an accurate 
evaluation of the hyperfine matrix elements is still unavailable. Overall, the fermion fit 
is the better of the two. This is more obvious at high fields when the anisotropy Ai — A2 
is high (ie. when the field is along the c-axis). For both models, the fit to the H\\b 
data becomes progressively worse as the field is increased. Fitting to the lower field data 
seems to give better overall agreement than fitting to the higher field results. This is not 
the case for the h\\c data (at least with the fermions). Since the field in the experiment 
was not actually placed along the chain c-axis, we might expect even worse agreement 
between this set of data and our calculations! In fact, as mentioned before, we expect a 
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very weak field dependence for tlie h\\b data wliicli would result from being close to U{1) 
symmetry. This was the universal result of Chapter 3. 

As is evident from the figures, the slope of the h\\b data and the calculated results 
agree. This implies that the relaxation is largely mediated through thermal bosons and 
that the calculation is off by a T-independent multiplicative factor. For small anisotropy, 
h ^ E, this effect cannot come from the matrix elements for the transition or the 
density of states. We believe that we have taken account of the obvious mechanisms for 
relaxation. Terms coming from the structural properties of NENP into the Hamiltonian 
(as discussed in the last section) are too small to be responsible for such a large increase 
in the relaxation at the mid-field range. Moreover, they would be expected to play 
a similar role when the field is placed along the c-axis. There are 16 protons in the 
chelate surrounding each Ni ion. Nuclear dipole-dipole interactions among them are 
energetically negligible, and thus could not be the cause for the increase in relaxation. 
It is certainly conceivable that the averaged hyperfine coupling is highly anisotropic, but 
it's hard to explain why there would an additional dependence on the magnitude of the 
field. Perhaps the discrepancy is due to reasons intrinsic to the experiment. 

Next we attempt to fit to the low field measurements taken for field along the 6-axis. 
We find that for fields less than 4 Tesla, it is not sufficient to consider the bulk theory 
alone. The relaxation rate decreases with increasing field in this regime (see Fig. 5.10| ). 
We can try to apply the impurity model to explain the data. Assuming the phase shift 
constants, A!|_, in Eqn. ( [4.10|) are 0(1), the impurity resonance width, P, derived in 



the last chapter can be graphed as in Fig. |5.11| . As is clear from the plot and Eqn. 
( [4.26| ), the impurity relaxation rate is essentially one delta-function peaked at D'/2 = h 
and another peaked at D' /2 + J' = h. This means that we expect two bumps in the 
relaxation rate due to impurity effects. The width of the bumps should correspond to 
the width of distribution of impurity couplings. The problem arises when we see that the 
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Figure 5.10: Relaxation rate for field along the 6-axis. 
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Figure 5.11: T{T) — the width of the impurity resonance 

temperature dependence of the low field data is roughly exponential: ~ e""^'/"^, where 
the impurity coupling J' is about 4.7K. Furthermore, the sharp decrease from zero field 
relaxation suggests D' = 0. By analyzing Eqn. ( [4.26| ) we see that the second bump 



should have little temperature dependence. This means that assuming the first bump 
sits near h = 0, the second must be larger and separated by about 3.5T. This is clearly 
not the case. Indeed, we would need a complicated distribution of couplings, J' and 
D', to get a proper fit. Adding an E type anisotropy will not change these conclusions. 
We thus do not have a satisfactory explanation for the low field behaviour. One should 
take notice, however, that the data was taken for a field along the 6-axis, where other 
problems were present at mid-field. 

Finally, we would like to mention some recent NMR data collected on the 1-D S = 1 
spin chain AgVP2SQ by Takigawa et. al. [Q. This material is highly one dimensional 
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with a large gap (A ~ 320K) and very nearly isotropic {5 ~ 4K). These characteristics 
make it ideal for analysis using our results. There are, however, some questions about the 
properties of the material which would have to be analyzed before an understanding of 
the NMR results is possible within the framework proposed here. The gap deduced from 
studies on the Vanadium atom (A ~ 410K) conflicts significantly with those performed on 
the phosphorus sites and with neutron scattering data. In addition, the material has very 
low symmetry (corresponding to the space group P2/a) and very little is known about 
the possible small E and D terms in the Hamiltonian and their corresponding symmetry. 
There is fair quahtative agreement between the ^^P NMR data and our theory, and it is 
possible to explain some of the discrepancies using a temperature dependent anisotropic 
gap structure, but we feel that not enough is yet understood about gross features of the 
material to justify such speculation at this time. 



Chapter 6 

Suggested Experiments and Curious Predictions 

We finish by pointing in this final chapter towards further experimental work which could 
serve to both better understand S — 1 IDHAF's as well as corroborate and clarify some 
of the issues raised in this thesis. 

6.1 Experimentally Testable Conflicts Between Models 

When discussing the matrix elements, < k,a\S\0)\q,b >, within the different models, 
we noticed that there were some discrepancies between predictions. We now examine 
this hoping to offer experiments that would resolve the issue in favour of one model or 

another. 

We start by discussing experiments on isotropic systems. In this case, the major 
differences between the predictions of the models concern large 9 transitions, where we 
recall from Chapter 2 that 

<k,a\r{0)\q,b>^ie'''^G(e) cosh(e) ^ -{ujk(^g - v^kq)/A'^ (6.1) 

This is especially dramatic in the case of backscattering. The problem with an experi- 
ment which probes large 6 transitions is that contributions from matrix elements of the 
staggered part of the spin may be large as well. This can be cured by looking for a low 
temperature experiment (T <S A), where the energy exchanged with the probe is small. 
As shown in the analysis of 1/Ti, the staggered contributions will be suppressed by a 
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double Boltzmann factor, e~^^/"^. A good candidate for such an experiment^ is elastic 
neutron scattering at zero or near zero magnetic field. The cross section is proportional 
to the spin correlation function; for elastic scattering, this is 

5(g,0) oc^l <n\ 5(0)|m> \^S{uJn - ujJS{Q - K + kje-'^"^/^ (6.2) 

n,m 

At sufficiently low temperatures, this expression is simpler than the analogous one for 
the relaxation rate thanks to the momentum conserving delta function. The energy 
conserving delta function ensures that only backscattering will contribute to the cross 
section. Using the results of Chapter 3 we easily integrate this to give 

5(Q,0) oc |G'(^)P^^e-"«/2/^ (6.3) 
vQ 

The NLa model gives 

tanh(^/2) 



e/2 



(6.4) 



64|l + (^/27r)2| 

At large Q this will behave as 1 / \og^ (vQ / A) . This is very different from the free boson 
prediction of G{6) = 1, and from the free fermion prediction of G{6) — > A'^/{vQ)^ for 
large Q. We need to qualify what we mean by 'large' Q. As discussed in Chapter 
1, the field theoretic models introduced are expected to be accurate only for Q near 
zero and vr. If we want to explore the two magnon nature of the structure function, 
we must be near Q ~ 0. What we mean by 'large' momentum elastic scattering is 
the investigation of the structure function near the border region where the field theories 



begin to diverge from numerical simulations [16 ; a region which satisfies all the criterion is 



.27r < Q < Att. This corresponds to energies three to six times that of the gap. We expect 
that the differences between the models should be discernible in this range. The reason we 
suggest the experiment be done at zero or nearly zero magnetic field is to ensure that only 



^T^ ^ relaxation is not an appropriate tool since the transitions arc dominated by small momentum 
transitions 
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backscattering transitions contribute. For nonzero field, interbranch transitions can occur 
at large momentum which will not necessarily select only backscattering events. This will 
not serve to make the interpretation transparent. The condition for backscattering even 
in the presence of a magnetic field is 



In the case of axial symmetry, we can suggest the same technique to investigate the 
difference between the zero field predictions of the boson and fermion models. Regardless 
of the size of D, if one only considers the cross section for scattering with Qv > J, 
then the fermion model predicts a result that vanishes as ~ ^ ^2^^ while the boson 
model prediction only involves the exponential factor. The same comments apply to the 
case where axial symmetry is broken as well. This is no surprise since at large enough 
momentum, 0(3) symmetry is effectively restored. 

Elastic neutron scattering is a good probe for the matrix elements involving large 
momentum and small energy exchange. Other techniques which explore the opposite 
regime are electron spin resonance (ESR) or far infrared absorption experiments. In both, 
one subjects the magnetic system to an external source of electromagnetic radiation (the 
microwave frequency value of the radiation depends on the transitions one is interested 
in investigating). The RF field couples to the spins in the same way that a magnetic field 
does, assuming that the electric dipole moment of the electrons on the magnetic ion is 
much smaller than the effective spins. ^ The interaction Hamiltonian is therefore 



Since the coupling is to the total spin of the system, the resonant transitions implied 

by Fermi's Golden Rule will involve energy uj and zero momentum exchange. At low 

rigorous treatment would try to treat the coupling to the electric dipole moment; this can be done 
within the spin manifold using the Wigner-Eckart theorem. We will not bother with such a treatment 
here, but we note that it may be crucial in understanding some experiments on NENP |32] 



Q > h/v 



(6.5) 



— » X — ^ — * 

Hj = H^Y^ ■ 2^ Si cos{ut) 



(6.6) 
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temperatures, the power absorbed when a uniform field is apphed to the system will be 
I{u) oc ^ I < n\Si^Q\m > \'^S{uJn{h) - uj^{h) - uj)6{kn - k^)e-'''^'^ 

n,m 



< a;Ku;l{h)\Sl^,\b;Ku;l{h) > \'e 



(6.7) 



dk 

where a and b denote one magnon states and k satisfies, uj^{h) — uiKh) = lo. Since 
the density of states factor in the above is divergent for A; = 0, it stands that /(a;) will 
have a peak at the value of h for which A"(/i) — t^iK) — uj. (The divergence will be 
cured by higher dimensional effects as discussed previously.) In a typical experiment, 
one judiciously chooses the RF frequency, a;, to be in the vicinity of desirable transitions, 
and the uniform field is then tuned to the peak in the absorption power. This is much 
easier to do than to fine tune the RF field. 

— * 

Let us now relate the ESR matrix elements to la,b{^i 0), calculated in Chapter 2. 

< a;0|5'^^ol^;0 >= j dx < a;Q\S\x)\b;Q >= J dx < a;0\e'^''S'{0)e-'^''\b;Q > 

= Jdx<a; 0\S'(0)\b; >^ Llib(0, 0) (6.8) 

Interesting confiicts between the models can be seen when there is some kind of, preferably 
large, anisotropy. For example, considering axial symmetry with a large D anisotropy. 



IC(0'0)P 



i(A3/Ax + Ax/A3 + 2) Bosons ^ ^ 

^ (6.9) 

2 Fermions 



The maximal difference corresponds to A^/ A± ~ 2 which leads to a discrepancy of about 
13% between the models. The closer the two branches lie, the better the agreement be- 
tween the models. This suggests the following experiment on highly anisotropic materials 
(NENP being a prime candidate). One chooses two RF frequencies. The first should cor- 
respond to the large interbranch gap, D, and the peak absorption ought to be measured 
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with a low field placed along the direction of the D anisotropy. The second RF frequency 
should be 0{E) if the material breaks axial symmetry, or 0{h) if the material is axially 
symmetric. This should then be used to measure the absorbed ESR power with the posi- 
tions of the uniform and RF fields exchanged. This second transition will involve matrix 
elements which will be gap independent in both models. The matrix elements from the 
first transition can be extracted and compared to that of the first. If the boson model 
is a better description even at these low fields, then the two matrix elements should be 
identical. 

One may argue that it is redundant to make both measurements since, if the gaps 
are known, Eqn. ( |6.7| ) should give the correct description. The problem lies in cutting 
off the infrared divergence at the absorption peak. This will introduce an unknown 
proportionality constant. This divides out when comparing the two measurements. The 
ratio of the two measurements would be 



3 



2 



/o+ (0,0) (6.10) 



/(cus) 1 - 

where we assume a small field, h -C Aj^, and small E D. 

To end this section, while on the subject of ESR experiments, we would like to propose 
additional experiments to test the impurity model presented in Chapter 4. ESR is ideal 
for such tests. Used in conjunction with Tf^ measurements on a given sample, it would 
be possible to characterize the couplings J' and D' of the end spins . 

6.2 Measuring Small Anisotropies 

Recall that we expect a peak in Tf^ whenever two branches cross. Experiments on 
Haldane Gap materials have yet to look for these. The sharpness of this peak depends 
on the interchain couplings which cut off the diverging integral in the calculation of 
the relaxation rate. Often, this will be broad because intrabranch transitions will share 
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the same cutoff (ie. wlien Jj > lun)- However, tlie bump sliould be experimentally 
observable. We propose that information about the anisotropy tensor can be extracted 
from this phenomenon. Essentially, one looks for the lowest field at which this bump 
occurs. This would give the direction of the lowest branch and the size of the anisotropy. 
We now explain this further. 

We assume the material in question has a well resolved D anisotropy and a seem- 
ingly degenerate doublet unresolved by other experimental techniques, such as neutron 
scattering. One begins by placing a uniform magnetic field in the plane perpendicular to 
the axial direction (ie. somewhere in the xy-plane). The magnitude of the field should 
be /i^ > 5(A3 — Apgrp), where 5 is the uncertainty in resolving the doublet. One then 
proceeds to measure Tf^ for different angles in the xy-plane spanning a region of at 
most 180°. If there is an E type anisotropy, one ought to see some structure to the data 
as a function of angle. Moreover, if there is such structure, we expect a bump at the 
angle where the branches cross. Once this angle is found, the experiment is repeated for 
somewhat lower field. The angle where the new bump should be seen would be greater 
than the old. There is actually enough information in these two measurements already 
to determine the anisotropy tensor. The dispersion relations are a function of the angle 
of the field (relative to some axis), the field magnitude and the gaps. The only unknowns 
are the absolute angle (or location of the axes of the anisotropy tensor) and the difference 
in gaps, I Aj; — Ay|. The two measurements could be used to solve for these two unknowns. 

In principle, one could also continue lowering the field and looking for the bump angle 
until it's clear that signal is being lost when the field is reduced further. At this point, 
one has located the minimum crossing field which must lie along the direction associated 
with the middle gap. This field also gives the anisotropy: |Aj. — Aj^| = /i^/jA^) — A|. 

It would be interesting to perform such an experiment on NENP. Presumably, one 
would find two angles corresponding to the two inequivalent chains in each unit cell. 
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Moreover, one would be able to verify the claims made in the last chapter regarding the 
positions of the local anisotropy tensor in NENP. 

6.3 ESR for NENP 

In the last chapter we noted that NENP has two inequivalent chains per unit cell. Fur- 
thermore, their local anisotropy tensor was argued to have symmetry axes which did not 
correspond to the crystal axes. These facts have important ramifications for ESR ex- 



periments on NENP. Figure |6.1| shows the dispersions for chains 1 and 2 (bold and light 
lines, respectively) when the field is 7r/3 from the crystallographic c axis in the ac-plane. 
This is an example of how transitions at two field strengths ought to be possible in the 
ESR experiment. 

Figure |6.2| shows the resonance field versus orientation of field in the crystallographic 
ac-plane. The lower branch denotes transitions in chain 1 while the upper branch cor- 
responds to transitions in chain 2. The transitions were calculated at .19 meV. corre- 



sponding to 47 GHz. In addition the experimental results of Date and Kindo are 
represented by the circles. One immediately sees that the data does not compare well 
with the predictions based on the models we've used so far, for instead of following one of 
the branches, the experimental results lie between them. Furthermore, it seems unlikely 
that perturbations will cause such a significant shift in the resonance field. One sees that 
the discrepancy is ~ ±1 Tesla. One possible explanation is that since the ESR signal in 
was also ~ ±1 Tesla in width and symmetric (in confiict with the predictions of ||47||), 



the signal from the resonances in both chains was somehow smeared and interpreted as 
one single peak. Seen that way the model predictions are in good agreement except for 
the large field regime. One also has to keep in mind that the high-field boson dispersions 
are not accurate and therefore the predictions at larger angles could easily be .5 Tesla 
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Figure 6.1: Dispersions for the two chain conformations and sample resonant transitions 
for a uniform field placed 60° from the crystallographic c-axis. 
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(or more) off the mark. We propose that further ESR experiments be done on NENP 
which specifically look for the double resonance predicted here. 

To end this discussion, we'd like to elaborate on a previously made statement re- 
garding the assignment of masses to the local Ni symmetry axes. It's easy to see that 
switching the masses around is tantamount to a 11/2 shift in Figure |5]^ (the fact that 
the gyromagnetic constants are not the same in orthogonal directions will not change 
the ESR resonance graph much since the ratio of the gyromagnetic constants is 0.98). 
Redrawing Figure |6.2| with this geometry misses the experimental results by 4 Tesla at 
and 90 degrees, where the two chain resonances coincide. This determines the proper 
labeling of the local symmetry axes. 



Chapter 7 



Concluding Remarks 



With the increasing theoretical interest in low dimensional systems, there has been a 
proportionate increase in the number of both realizable physical systems and experiments. 
This work offers a comprehensive analysis of NMR relaxation in Haldane gap materials, 
taking account of anisotropy and other material properties. As well, our analysis has led 
to predictions pertaining to other types of experiments. It is hoped that our efforts will 
aid in both extending and clarifying existing knowledge of the subject. 
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